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ABSTRACT 


Experiments  were  carried  out  using  an  existing  tube 
alone,  and  with  the  tube  connected  to  a  two-dimensional 
wedge  nozzle.  The  range  of  maximum  duration  of  steady 
reflected  pressure  from  3.5  to  5  milliseconds  was  achieved 
through  tailored  operation  for  incident  shock  strengths  of 
3.4  and  2.0,  using  pure  Helium  and  a  70%  Helium/30%  Nitrogen 
mixture  as  the  driver  gas  respectfully. 

Spark  and  continuous  light  shadowgraph  techniques  were 
attempted  using  an  optical  window  at  the  Mach  4.3  location. 
Results  demonstrated  that  the  short  duration  flow  phenomena 
in  a  shock  tunnel  can  be  recorded  successfully  using 
existing  equipment.  Calculations  showed  that  the  addition 
of  a  Mach  10  nozzle  and  15  m^  (6'  diameter  x  15'  long)  dump 
chamber  would  provide  a  useful  hypersonic  facility  for 
instruction  and  research. 
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I. 


INTRODUCTION 


The  simulation  of  flight  conditions  using  ground  based 
facilities  has  been  an  essential  part  of  all  aircraft 
development.  Figure  1  illustrates  the  various  types  of 
facilities  which  have  evolved  to  simulate  flight  at 
progressively  higher  speeds  and  altitudes.  Increasing 
energy  requirements  led  to  blow-down  wind  tunnels.  The  need 
to  avoid  condensation  in  nozzle  expansions  required  heaters. 
The  very  high  stagnation  enthalpies  associated  with 
hypersonic  flight  led  to  very  short  duration,  intermittent 
techniques.  Each  of  the  types  shown  in  Figure  1  was  in 
operation  more  than  20  years  ago.  A  renewed  interest  in 
hypersonic  simulation  has  been  generated  by  the  National 
Aerospace  Plane (NASP)  program. 

A  recent  review  of  facility  requirements  for  hypersonic 
propulsion  system  testing  by  Dunn  and  Lori  [Ref.  1] 
concludes  that  a  single  ground  based  facility  will  not  be 
capable  of  providing  the  necessary  testing  to  cover  the  NASP 
flight  envelope.  The  planned  use  of  Computational  Fluid 
Dynamic  (CFD)  codes  to  bridge  the  gap  left  by  inadequate 
ground  based  facilities  will  make  CFD  codes  an  essential 
part  of  the  design  process  but,  before  confidence  can  be 
placed  in  these  codes,  fundamental  experiments  are  needed  to 
validate  the  models  incorporated  for  real  gas  behavior. 
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Figure  1.  Ground  Based  Testing  Facilities  [Ref.  1] 
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Recognizing  the  need  for  CFD  validation,  the  feasibility 
of  establishing  a  shock  tunnel  capability  by  modifying  the 
existing  Naval  Postgraduate  School  (NPS)  three  inch  diameter 
shock  tube,  was  examined  in  the  present  study.  The  time 
duration  and  magnitude  of  constant  pressure  available  in  the 
reflected  region  of  the  shock  tube  were  evaluated.  The 
concept  of  tailoring  was  examined  using  two  types  of  driver 
gases,  pure  Helium  (He)  and  a  mixture  of  70%  Helium  and  30% 
Nitrogen  (N2)  ,  driven  into  Nitrogen  as  a  test  gas.  An 
exploratory  shock-tunnel  nozzle  experiment  was  conducted, 
and  the  capability  of  recording  short-duration  flow 
phenomena  using  available  visualization  techniques  was 
demonstrated . 

It  was  concluded  that  the  development  of  an  operational 
hypersonic  shock  tunnel  was  feasible  and  would  provide  the 
means  for  follow-on  efforts  in  the  area  of  CFD  code 
validation.  At  the  same  time,  such  a  facility  would  provide 
a  useful  illustration  of  hypersonic  flow  to  aeronautical  and 
space  applications  engineering  students. 

The  following  material  is  divided  into  three  sections; 

1)  a  review  of  shock  tube  and  shock  tunnel  ideas,  2)  a 
discussion  of  the  experimental  investigation  including  a 
description  of  the  facilities,  experimental  procedure  and 
results;  and  (3)  an  evaluation  of  nozzle  throat  and  dump 
chamber  sizing  requirements  for  a  hypersonic  tunnel. 
Instrumentation  modifcations  will  be  required  in  the 
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development  of  a  hypersonic  capability.  Appendix  A  reviews 
the  existing  facility  and  provides  recommendations  for 
upgrading  the  current  hardware  and  instrumentation. 
Simplified  flow  models  are  developed  in  Appendices  B  and  C 
to  enable  the  selection  of  nozzle  throat  diameter  and  dump 
chamber  volume  respectively.  Experimental  pressure  and  thin 
film  traces  obtained  in  the  experimental  program  are 
contained  in  Appendix  D  together  with  predicted  results 
obtained  using  a  computer  performance  code,  which  is  listed 
in  Appendix  E. 
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II.  THEORETICAL  CONSIDERATIONS — A  REVIEW 

A.  BASIC  SHOCK  TUBE 

A  shock  tube  is  an  economical  short  duration  test 
facility  used  in  laboratory  experiments  in  gas  dynamics, 
aerodynamics  and  physics.  The  principles  of  shock  tube 
operation  ,  which  are  well  documented  in  numerous  texts,  for 
example,  Gaydon  and  Hurle  [Ref.  2]  and  Owczarek  [Ref.  3], 
are  illustrated  in  Figures  2  and  3  [Ref.  4].  As  shown  in 
Figure  2,  the  device  consists  of  a  long  tube  divided  into  a 
high  pressure  and  a  low  pressure  section  by  a  diaphragm. 

When  the  diaphragm  separating  the  high  pressure  driver 
section,  region  4,  and  the  low  pressure  section,  region  1, 
is  ruptured,  a  series  of  compression  waves  coalesce  into  a 
single  shock  front  (S)  which  compresses  and  heats  the  region 
1  test  gas  to  region  2  conditions.  Simultaneously,  a  series 
of  rarefaction  (expansion)  waves  (R)  travel  into  the  high 
pressure  section,  region  4,  cooling  and  expanding  the  driver 
gas  to  region  3  conditions. 

Region  3  (the  expanded  driver  gas)  and  region  2  (the 
compressed  test  gas)  are  separated  by  a  contact  surface  (C) . 
The  gases  on  either  side  of  the  contact  surface  travel  with 
the  same  velocity,  and  are  at  the  same  pressure,  but  are 
different  in  temperature.  A  time  history  of  the  process 
(time  vs.  distance)  is  shown  in  Figure  3. 
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Figure  3.  Time  History  of  Events  after  Diaphragm 
is  Broken  [Ref.  4] 


What  is  shown  in  Figure  3  can  be  determined  accurately 
using  the  method  of  characteristics  and  shock  theory  (see 
[Ref.  3]).  As  the  shock  front  reflects  off  the  driven  tube 
endwall,  the  region  2  gas  is  effectively  shocked  a  second 
time,  creating  in  region  5  stagnant  gas  at  high  pressure  and 
temperature.  This  region  may  be  used  for  high  temperature 
gas  dynamic  studies  or,  when  coupled  to  a  nozzle,  made  to 
serve  as  a  reservoir  of  compressed  and  heated  gas  for  a 
hypersonic  wind  tunnel . 

B .  PERFORMANCE 

The  strength  of  the  generated  shock  (M^)  is  related  to 
the  pressure  ratio  (P^/PJ  across  the  diaphragm  and  the  speed 
of  sound  in  the  undisturbed  gases  through  the  equation  given 
by  Gaydon  [Ref.  2]. 
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2,iM3-(>,-1) 


(1) 


Figures  4  [Ref.  5]  and  5  [Ref.  6],  which  illustrate  Equation 
(1) ,  show  that  the  strongest  shocks  are  generated  by  using 
high  pressure  ratios  and  by  using  driver  gases  with  high 
speeds  of  sound  and  low  ratios  of  specific  heats. 

As  the  desired  shock  strength  increases,  the  required 
P^/Pi  increases  significantly.  To  achieve  stronger  shocks 
several  approaches  are  possible.  The  most  common  include 
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gure  4 . 


Pressure  Ratio  (P4/P1)  Required  to 
Generate  a  Given  Shock  Strength  Using 
Various  Gas  Mixtures  [Ref.  5] 


(b) 


Driver  Gas  Specific  Heat  of  1,67  and  a  Driven 
Gas  Specific  Heat  of  1.67 


Figure  5  (CONTINUED) 


the  use  of  an  area  reduction  between  the  driver  and  driven 
sections  [Ref.  5],  the  use  of  multiple  diaphragms  with  a 
buffer  gas  in  the  driver  [Ref.  5],  and  heating  the  driver 
gas  through  combustion  or  electrical  discharge  techniques 
[Ref.  6].  Figures  6  and  7  [Ref.  7]  show  the  relationship 
between  shock  strength  (M,)  and  reflected  region  5 
conditions  for  equilibrium  air. 

C.  SHOCK  TUNNELS 

The  reflected  region  5  pressure  and  temperature  may  be 
used  to  supply  a  supersonic  or  hypersonic  nozzle  if  the 
nozzle  throat  is  made  sufficiently  small  to  ensure  full 
reflection  of  the  incident  shock.  The  evolution  of  this 
arrangement  is  shown  schematically  in  Figure  8.  To  maximize 
test  time  and  minimize  the  nozzle  starting  time,  often  a 
diaphragm  separates  the  nozzle  and  dump  chamber  from  the 
driven  end  of  the  shock  tube.  The  dump  chamber  is  sized  and 
evacuated  to  a  level  which  ensures  that,  when  the  tube  is 
fired,  the  pressure  in  the  chamber  remains  below  the  static 
pressure  at  the  nozzle  exit  for  the  duration  of  steady  flow. 
[Ref  8:pp.  624-640] 

For  hypersonic  applications  a  conical  nozzle  is  often 
used.  Table  1  provides  calculated  test  conditions  for 
various  flow  expansions  using  ideal  gas  theory.  The 
computer  program  used  for  the  calculations  is  described  in 
Appendix  E.  Since,  in  an  isentropic  expansion,  static 
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Figure  6.  Reflected  Pressure  as  a  Function  of  Shock 
Strength  for  Eguilibrium  Air  [Ref.  7] 
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Figure  7.  Reflected  Temperature  as  a  Function  of  Shock 
Strength  for  Equilibrium  Air  [Ref.  7] 
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Figure  8.  Evolution  of  a  Shock  Tunnel 
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TABLE  1 


THEORETICAL  PREDICTION  OF  TEST  CONDITIONS  FOR  A  HYPERSONIC 
NOZZLE  (Me  =  10)  DRIVEN  FROM  THE  REFLECTED  REGION 

OF  A  SHOCK  TUBE 


Driver 

P5 

T5 

Pe 

Te 

^start 

Gas 

Ms 

(Atm) 

(°K) 

(mm  Hg) 

(°K) 

(msec) 

“2.4 

58 

973 

1.02 

50 

.041 

3.7 

32 

2000 

.55 

94 

.  029 

Helium 

4.2 

23 

2460 

.41 

116 

.  026 

4 . 5 

19 

2700 

.33 

131 

.024 

""  1.8 

23 

632 

.4 

30 

.  051 

2.4 

9.7 

971 

.  169 

50 

.  041 

Nitrogen 

2 . 6 

6.7 

1110 

.118 

52 

.038 

2.8 

5.3 

1200 

.  09 

56 

.  037 

temperature  decreases  as  the  Mach  number  (of  the  test 
section)  increases,  liquefaction  of  the  test  gas  becomes  an 
issue.  Since  nitrogen  closely  matches  the  aerodynamic 
properties  of  air,  but  liquefies  at  a  lower  temperature,  it 
is  preferred  for  hypersonic  facilities  in  which  simulated 
real  gas  effects  are  not  required.  This  is  illustrated  in 
Figure  9  which  also  shows  the  interesting  results  of  wind 
tunnel  experience  [Ref.  9].  In  practice,  condensation  is 
observed  to  occur  at  temperatures  which  are  lower  than  the 
saturation  temperature,  thereby  reducing  somewhat  the 
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Pressure,  mm  Hg 


Figure  9.  Experimentally  Determined  Liquif ication 
Properties  of  Air  and  Nitrogen  [Ref.  9] 


stagnation  temperature  requirements  to  reach  a  given  Mach 
nuiaber . 


Given  the  extremely  short  test  time  available  in  a  shock 
tube,  any  coupling  to  a  hypersonic  nozzle  must  be  done  in 
such  a  way  so  as  to  minimize  nozzle  starting  transients  and 
provide  the  maximum  available  steady  flow  in  the  test  sec¬ 
tion.  The  starting  time,  defined  as  the  time  to  approach 
steady  flow  through  a  hypersonic  conical  nozzle  (71  =  7/5)  , 
is  given  by  Smith  [Ref.  8]  as 
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A  reduction  in  starting  time  may  be  achieved  by  using  a 
large  nozzle  angle  (fi  =  nozzle  half  angle) ,  a  small  throat 
radius  (R^)  and  a  high  stagnation  temperature  in  the 
reflected  region.  [Ref.  8:pp.  625-640]  The  variation  in 
starting  time  with  increasing  incident  shock  strength  for  a 
conical  nozzle  with  fi  =  30  degrees  is  shown  in  Table  1. 

After  the  unsteady  starting  process,  the  useful  test 
flow  will  continue  as  long  as  a  sufficiently  high  pressure 
ratio  is  maintained  across  the  nozzle.  The  steady  flow 
process  will  be  a  function  of  only  the  upstream  reflected 
stagnation  conditions  as  long  as  the  dump  chamber  has 
sufficient  volume.  Ideally,  uniform  reservoir  stagnation 
conditions  will  end  when  the  reflected  region  conditions  are 
changed  by  a  wave  emanating  from  the  interaction  of  the 


18 


contact  surface  and  the  head  or  tail  of  the  reflected  expan¬ 
sion  wave  (Figure  2) .  The  useful  test  time  in  the  nozzle 
ends  when  the  reflected  pressure  drops  below  the  level 
rec[uired  to  maintain  the  static  pressure  at  the  nozzle  exit 
above  the  back  pressure  in  the  dump  chamber.  Thus  the 
useful  test  time  will  depend  on  the  nozzle  throat  size, 
which  governs  the  flow  rate  out  of  the  tube,  and  the  size  of 
the  dump  chamber,  which  determines  the  rate  at  which  the 
back  pressure  builds  up. 


D.  TAILORED  OPERATION 
1.  Tailoring 

The  fundamental  limitation  of  a  shock  tube  is  the 
available  test  time.  Considerable  effort  has  gone  into 
trying  to  extend  shock  tube  test  times  [Ref.  10].  Test 
times  can  be  lengthened  by  as  much  as  a  factor  of  ten  by 
"tailoring,”  which  means  adjusting  conditions  so  that  the 
ratio  of  the  speed  of  sound  and  specific  heat  behind  the 
reflected  shock  matches  those  in  the  region  behind  the 
contact  surface.  When  this  tailored  condition  is  present, 
the  reflected  shock  passes  through  the  contact  surface 
bringing  it  to  rest.  A  Mach  wave  is  reflected  forwards  off 
the  contact  surface  and  generates  no  additional  disturbances 
as  it  passes  through  region  5.  In  the  untailored  condition, 
the  wave  reflected  forward  from  the  contact  surface  is  not  a 
Mach  wave  but  a  rarefaction  wave  when  a,  <  aj,  and  a  shock 


wave  when  ai  >  a^.  The  reflected  waves  cause  conditions  in 
region  5  to  change.  Experimental  evidence  of  the  effect  of 
the  wave  interactions  obtained  by  Pennelegion  [Ref.  10]  are 
shown  in  Figure  10.  Pennelegion' s  research  fully  explored 
tailoring  analytically  and  provided  supporting  experimental 
evidence.  His  findings  and  experimental  techniques  were 
followed  in  the  experiments  undertaken  in  the  present  study. 
The  tailored  interface  technique  requires  that  the  following 
equation  be  satisfied  [Ref.lOrp.  4]; 
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A  unique  tailored  condition  (shock  Mach  number) ,  exists  for 
every  combination  of  driver  and  driven  gas.  The  duration  of 
constant  pressure  in  region  5  lasts  until  the  rarefaction 
wave  reflects  from  the  driver  endplate  and  returns  to  meet 
the  contact  surface.  For  low  incident  shock  strengths  it  is 
the  head  of  the  rarefaction  wave,  and  for  high  incident 
shock  strengths  it  is  the  tail  of  the  rarefaction  wave  which 
terminates  the  constant  pressure. 

The  tailored  interface  technique  for  producing 
extended  periods  of  constant  reflected  pressure  and 
temperature  was  demonstrated  experimentally  for  Mach  numbers 
below  3.4.  When  the  technique  was  applied  to  shocks  of  Mach 
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Figure  10. 


Experimental  Observation  of  the  Three  Types 
of  Wave  Interactions  in  the  Reflected  Region 
of  a  Shock  tube  [Ref.  10] 


number  greater  than  3.8,  there  was  a  significant  drop  in 
test  time  due  to  boundary  layer  interaction.  In  the  under¬ 
tailored  condition  the  reflected  conditions  showed  an 
initial  fall  in  temperature  and  pressure,  followed  by  a  slow 
rise  in  pressure  and  continued  drop  in  temperature.  The 
overtailored  condition  revealed  a  pressure  rise  from  the 
reflection  from  the  contact  surface  followed  by  a  rapid  drop 
in  temperature.  In  all  cases  with  strong  incident  shock 
waves,  the  real  gas  effects  caused  significant  variation  in 
the  reflected  region  conditions.  [Ref.  10] 

2 .  Mixing  of  Driver  Gases  THe  and 

The  application  of  tailoring  to  increase  test  time 
in  a  hypersonic  shock  tunnel  allows  more  accurate 
measurements  of  force  and  pressure  to  be  made  on  test 
models.  As  desired  test  section  Mach  number  increases  there 
is  a  significant  decrease  in  working  section  static 
pressure.  It  is  therefore  desirable  to  have  reflected 
region  pressures  as  steady  and  as  high  as  possible.  Pressure 
recovery,  defined  as  the  ratio  of  the  initial  step  in  the 
reflected  pressure  (P5)  to  driven  pressure  (PJ  ,  may  be 
increased  by  increasing  the  speed  of  sound  of  the  driver  if 
some  test  time  is  sacrificed.  The  combination  of  a  driver 
of  Helium  (He)  and  Nitrogen  (N2)  can  provide  an  improvement 
in  test  time  and  pressure  recovery  over  a  driver  of  N2 


alone.  Experimental  evidence  of  this  is  shown  in  Figures  11 
and  12  from  Pennelegion  [Ref.  10].  From  the  same  study. 


I  2  1  4  Mi 

Figure  11.  Available  Steady  Pressure  Time  [Ref.  10] 
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Figure  12.  Experimental  Pressure  Recovery  for  Different 
Driver  Gases  [Ref.  10] 
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Figure  13  and  Table  2  show  tailoring  using  various  mixtures 
of  He  and  N2  to  generate  shocks  in  the  range  1  <  M,  <  3.4. 

The  tailoring  shock  Mach  number  may  be  selected  to  provide 
the  required  stagnation  temperature  to  prevent  liquefaction 
and  maximize  test  time  for  hypers  ..nic  nozzles  with  various 
Mach  numbers.  Pennelgion  [Ref.  10],  provided  the  method  for 
calculating  tailored  conditions  for  He,  N2  mixtures. 
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Figure  13.  Variation  of  Tailored  Mach  Number  with 
Helium  Fraction  [Ref.  10] 
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TABLE  2 


EXPERIMENTALLY  TAILORED  MACH  NUMBERS  FOR  DRIVER 
GAS  MIXTURES  [REF.  10] 


He 

Fraction 

Mte 

Pa/ Pi 

Time  Steady 
Pressure  (msec) 

T5 

(°K) 

0.7 

2 . 0 

24 

6.5 

714 

00 

• 

o 

2.7 

60 

6.0 

1096 

0.9 

2.8 

70 

4.7 

1159 

1.0 

3.4 

124 

4 . 3 

1553 
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III.  SHOCK  TUBE 


In  this  section  the  existing  shock  tube  and  instrumenta¬ 
tion,  and  the  modifications  made  during  the  course  of  the 
present  experiments,  are  described. 

A.  STRUCTURE 

The  three  inch  diameter  shock  tube  at  the  Naval 
Postgraduate  School  Aeropropulsion  Laboratories  is  shown  in 
Figure  14.  A  detailed  description  is  given  in  Penaranda 
[Ref.  4].  The  tube,  constructed  of  highly  polished  303 
stainless  steel  consists  of  a  driver  section  7.75  feet  long, 
shown  in  detail  in  Figure  15,  and  a  driven  section  18.9  feet 
long. 

B.  DIAPHRAGM  CHAMBER 

The  tube  incorporates  a  double  diaphragm  arrangement  to 
maintain  precise  and  repeatable  diaphragm  bursting  pressure 
ratios  (pypj  .  The  tube  currently  is  fitted  with  diaphragms 
designed  to  burst  using  a  constant  driver  pressure  of  600 
psi.  Scored  brass  diaphragms  designed  to  rupture  at  400  psi 
are  inserted.  Variation  in  shock  speed  is  obtained  through 
adjustment  of  the  downstream  pressure  (Pj)  .  To  obtain  the 
desired  pressure  ratio,  the  2.5  inch  long  diaphragm  chamber 
(Figures  15  and  16)  and  the  driver  chamber  are  simul¬ 
taneously  filled  to  one  half  of  the  desired  P4  pressure  (300 
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Figure  16.  Dual  Diaphragm  Chamber  (Unbolted) 

psi) .  A  diverting  valve  is  closed  to  isolate  the  diaphragm 
chamber  at  300  psi.  The  driver  section  pressure  is  then 
raised  to  the  desired  P<,  pressure  (600  psi).  To  fire  the 
tube,  a  diverting  valve  is  opened  and  some  of  the  driver  gas 
at  600  psi  is  dumped  into  the  diaphragm  chamber  causing  the 
first  diaphragm  to  break.  The  second  diaphragm  breaks  im¬ 
mediately,  as  the  pressure  on  the  driven  face  of  the 
diaphragm  rapidly  falls. 
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The  dual  diaphragm  arrangement  minimizes  bursting 
pressure  fluctuations  due  to  inconsistencies  in  the  material 
of  the  diaphragms.  By  applying  a  pressure  of  200  psi  above 
the  design  bursting  pressure  of  the  first  diaphragm,  rupture 
of  the  two  diaphragms  is  ensured  with  a  precise  value  of  P^. 
In  shock  tubes  using  a  single  diaphragm  arrangement,  precise 
bursting  pressures  cannot  be  repeated  unless  the  complexity 
of  a  diaphragm  puncturing  device  is  used.  Without  a 
puncture  device,  the  bursting  pressure  fluctuations  remain  a 
function  of  the  uniformity  of  the  diaphragm  material. 

C.  NOZZLE  ASSEMBLY 

An  existing  two  dimensional  stainless  steel  nozzle 
designed  for  an  exit  Mach  number  of  5.72  (A/A*  =  43.7)  and  a 
rectangular  dump  chamber  to  which  it  bolted  were  available 
for  the  present  study.  The  assembly  is  shown  in  Figure  17. 
For  the  present  experiments  two  1”  diameter  window  ports 
were  installed  to  allow  optical  measurement  techniques  to  be 
exercised,  and  a  model  support  was  constructed  to  position  a 
3/8"  axisymmetric  model  on  the  nozzle  centerline  in  line 
with  the  window,  where  the  Mach  number  was  about  4.3.  The 
final  assembly  is  shown  in  Figure  18.  The  dimensions  of  the 
nozzle  are  given  in  Figure  19  and  the  design  of  the  window 
installation  is  shown  in  Figure  20.  Views  of  the  nozzle  and 
model  support  are  shown  in  Figure  21.  Since  the  nozzle  and 
dump  tank  existed  (from  a  previous  gas  dynamic  laser 
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iraensional  Nozzle  and  Dump  Chamber 


Figure  18. 


Assembled  Dump  Chamber,  Model  Support 
and  Nozzle 


Existing  Nozzle 


Modified  Nozzle 
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investigation) ,  location  and  sizing  of  the  window  as  well  as 
the  itiaxiitium  flow  duration  were  limited  by  the  existing 
hardware . 

D.  INSTRUMENTATION 
1 .  Shock  Speed 

Shock  speed  was  determined  by  measuring  the  time  for 
the  shock  to  travel  between  thin  film  gauges  located  at 
positions  A  and  B  and  from  B  to  the  pressure  transducer  at 
C  (Figure  14) .  The  locally  manufactured  thin  film  gauges 
shown  in  Figure  22,  consist  of  a  0.25”  diameter  pyrex  rod 
coated  with  a  1/32"  wide  thin  film  of  platinum  giving  a  cold 
resistance  of  100  ohms.  [Ref.  4:p.  23]  As  the  region  of 
hot  gas  immediately  behind  the  shock  front  passes  over  the 
platinum,  its  temperature  and  resistance  change  and  a 
voltage  pulse  is  generated.  In  the  present  work  the  signals 
from  the  thin  films  were  amplified  and  recorded  for 
analysis.  The  thin  films  at  location  A  also  served  as  a 
trigger  for  the  digital  counter  and  shadowgraph  spark 
source. 

Shock  speed  was  determined  in  different  ways.  The 
pre-existing  techniques  involved  using  oscilloscope  records 
of  thin  film  and  pressure  transducer  outputs,  and  also  the 
output  of  a  digital  counter  triggered  on  and  off  by  the  thin 
films  at  A  to  B  respectively.  In  the  present  work,  in  runs 
1  through  49  shock  speeds  were  determined  from  a  Polaroid 
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photograph  of  the  oscilloscope,  such  as  is  shown  in  Figure 
23,  and  from  the  digital  counter.  Both  techniques  were 
identical  to  those  used  in  the  late  1960 's. 

In  runs  50-85,  the  outputs  of  the  three  traces  were 
recorded  digitally  on  three  of  the  available  four  channels 
of  a  Scientific  Alanta  SD380  Spectrum  Analyzer.  Data  were 
recorded  at  .019  and  .039  msec  intervals  depending  upon  the 
time  scale  selected.  An  example  is  shown  in  Figure  24.  The 
shock  speeds  were  determined  during  postprocessing  of  the 
data . 

2  .  Driver  Ttibe  Pressure 

A  dial  gauge  with  10  psi  pressure  increments  was 
used  to  determine  when  the  constant  600  psi  driver  pressure 
was  reached.  A  dead  weight  tester  was  used  to  calibrate  the 
gauge . 

3 .  Driven  Tube  Pressure 

Two  calibrated  Wallace  and  Tiernan  Bourdon  tube 
gauges  were  used  to  determine  the  driven  tube  pressures.  In 
the  range  from  0  to  50mm  I!g,  a  gauge  with  increments  of  .2 
mm  Hg  was  used.  For  higher  driven  pressures,  a  second  gauge 
with  a  range  from  0  to  35  psi  in  increments  of  .05  psi  was 
used . 

4 .  Driven  Tube  Temperature 

A  Type  K  (chromel-alumel)  thermocouple,  connected  to 
a  Daigger  Scientific  Digital  Panel  Meter  using  matched 
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extension  wires,  was  mounted  in  the  driven  end  of  the  tube. 
The  digital  readout  was  recorded  prior  to  every  run. 

5 .  Reflected  Pressure  and  Test  Duration 

A  piezoelectric  pressure  transducer,  Endevco  model 

2501,  a  charge  amplifier,  Endevco  model  2700,  and  a 
Tektronics  551  oscilloscope  were  initially  used  to  record 
the  transient  region  5  pressure.  The  lack  of  resolution  and 
strange  qualitative  behavior  observed  in  the  results  (Figure 
23) ,  led  to  the  mounting  of  an  improved  transducer. 

A  PCB  Piezontonics,  Inc.  Model  124A  pressure  trans¬ 
ducer  connected  to  a  model  504A  dual  gain  charge  amplifier 
were  used  in  runs  50  through  75.  The  system  was  calibrated 
using  a  standard  static  calibration  feature  described  in  the 
manufacturer's  manuals  [Refs.  11,12].  The  transducer 
provided  a  response  of  240  psi/volts.  Experimental  traces 
obtained  are  given  in  Appendix  D. 

6 .  Oscilloscopes 

Two  Tektronics  oscilloscopes,  models  549  and  551, 
fitted  with  Polaroid  cameras  were  used  for  runs  1  through 
49.  The  capability  of  the  scopes  are  fully  described  in  the 
associated  manuals  [Refs.  13,14].  Time  measurements  to 
determine  shock  speed,  and  voltage  levels  to  determine 
pressure,  were  determined  by  overlaying  a  scale  on  the 
Polaroid  photographs. 
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7 .  Signal  Analyzer 


A  Scientific  Alanta  SD380  four  channel  signal 
analyzer  (Figure  25)  equipped  with  a  hard  disk  and  plotter 
was  used  to  record  and  process  data  from  runs  50  through  85. 
The  independent  channel  trigger  delays  (±) ,  adjustable 
scales  and  post-processing  capabilities  of  the  digitized 
data  (in  contrast  to  photographing  oscilloscopes) ,  is 
evidence  of  the  advances  in  measurement  capabilities  over 
the  past  20  years.  A  full  description  of  capabilities  of 
the  analyzer  is  given  in  the  instrument  manual.  [Ref.  15] 

E.  OPTICAL  MEASUREMENT  SYSTEMS 

A  simple  shadowgraph  arrangement  was  generated  from  the 
components  of  a  6"  schlieren  system.  The  optical  layout  is 
shown  in  Figure  26  [Ref.  16].  The  spark  source  was  used 
with  Polaroid  film  in  a  camera  back  attached  to  the  nozzle 
window.  The  continuous  source  was  used  with  a  Hycam  movie 
camera . 

1 .  Light  Sources 

a.  High  Pressure  Mercury  Arc  Lamp 

A  high  pressure  mercury  arc  lamp  (Per  Type  112- 
5145)  constant  light  source  was  used  with  a  six  inch 
achromatic  convex  lens  with  a  focal  length  of  60  inches.  The 
arrangement,  mounted  on  an  adjustable  bench  (Figure  27) , 
provided  a  parallel  light  beam  through  the  nozzle  windows. 
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A 


Figure  25.  Scientific  Alanta  SD380  Four 
Channel  Signal  Analyzer  with 
Hard  Dis.’:  Installed 
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6"  CONOENSTtlG  I.ENS 


-  SPARK  LIGHT  SOURCE 

-  MERCURY  ARC  LAMP 


Figure  26.  Optical  Layout  Sketch 
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Figure  27.  Constant  Light  Source  and  Lens  Mounted 
on  an  Adjustable  Optical  Bench 
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The  beam  was  adjusted  by  the  procedures  described  in  the 
associated  manual  [Ref.  16]. 

b.  NASA  229-D  Spark  Light  Source 

A  NASA  229-D  spark  light  source  (Figure  28) 
provided  a  nanosecond  spark  for  shadowgraph  pictures  taken 
with  Polaroid  film.  The  spark  was  triggered  from  a  thin 
film  connected  through  the  delay  and  voltage  circuits  shown 
in  Figure  29.  The  timing  for  the  spark  source  was 
determined  experimentally. 

2 .  Camera  Systems 

a.  Hycam  Movie  Camera 

A  Hycam  movie  camera  [Ref.  17]  with  a  framing 
rate  on  the  order  of  6000  pictures  per  second  at  runout  of  a 
400  ft  roll  of  film,  was  used.  The  camera  was  equipped  with 
a  Nikon  55mm  lens  and  ASA  250  black  and  white  film.  A  time 
reference  was  provided  on  the  edge  of  the  film  by  a  strobe 
light  pulsing  at  1000  pulses  per  second.  The  camera  was 
triggered  manually.  The  camera  was  mounted  on  a  second 
adjustable  bench  as  shown  in  Figure  30. 

b.  Spark  Shadowgraphs 

Spark  shadowgraphs  were  taken  using  a  Polaroid 
Land  Film  Holder  with  Type  52  ASA  400  Black  and  White  film, 
held  directly  behind  the  optical  windows.  To  eliminate 
stray  light,  the  mounting  arrangement  was  contained  in  a 
cardboard  box,  and  tests  were  conducted  after  dark. 
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Figure  29.  Spark  Timing  Delay  Circuit  and  High 
Voltage  Pulse  Generator 


Figure  30.  Hycam  Movie  Camera  with  Separate 
Pulse  Generator 
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IV.  EXPERIMENTAL  PROCEDURE 

A.  PREPARATION 

Prior  to  testing,  with  the  driven  endplate  removed  and 
the  diaphragm  chamber  opened,  a  cloth  lightly  dampened  with 
Acetone  was  pulled  through  the  driver  and  driven  sections. 
0-ring  grooves  of  the  opened  sections  where  thoroughly 
cleaned  and  a  thin  coating  of  vacuum  grease  was  applied  to 
the  rubber  0-rings.  Diaphragms  were  wiped  clean  of  visible 
dirt  and  the  system  reassembled.  The  cleaning  procedure  was 
repeated  after  approximately  15  runs.  Using  a  roughing 
pump,  a  clean-tube  vacuum  pressure  of  3  mm  Hg  was  routinely 
achieved. 

B.  LOADING  THE  DRIVEN  GAS 

Prior  to  loading  to  the  desired  test  pressure  (Pj)  ,  the 
driven  section  was  purged.  To  purge,  the  tube  was  evacuated 
to  3  mm  Hg  and  then  isolated  from  the  roughing  pump. 

Nitrogen  was  introduced  to  a  pressure  1  to  4  psi.  The  tube 
was  then  evacuated  a  second  time  and  repressurized  to  the 
desired  test  condition. 

C.  LOADING  THE  DRIVER  GAS  AND  FIRING 

The  driver  assembly  is  designed  to  operate  at  600  psi 
with  single  or  dual  gas  combinations.  To  purge  the  driver 
section  and  diaphragm  chamber,  the  system  was  pressurized  to 


approximately  30  psi  and  then  vented  to  atmospheric  pres¬ 
sure.  The  system  was  then  pressurized  to  600  psi  with  the 
diverter  valve  to  the  diaphragm  chamber  being  closed  at  300 
psi.  When  the  70%  He  and  30%  N2  driver  mixture  was  used, 
the  system  was  pressurized  to  420  psi  with  He  and  then  to 
600  psi  with  N2. 

The  tube  was  fired  by  opening  the  diverter  valve  to  the 
diaphragm  chamber.  Some  seepage  from  the  driver  section  was 
evident.  To  minimize  this  effect  the  tube  was  fired 
immediately  after  reaching  the  600  psi  level.  Replacement 
of  0-rings  would  minimize  this  problem.  Seepage  into  the 
driven  section  was  not  evident. 
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V.  EXPERIMENTAL  RESULTS  AND  DISCUSSION 


A  series  of  experiments  was  conducted  to  gain  experience 
with  the  shock  tube.  The  study  consisted  of  approximately 
85  shots  involving  the  use  of  He  and  a  70%  He/30%  N2  mixture 
as  driver  gases,  driven  into  a  N2  test  gas.  Studies  with 
the  shock  tube  assembled  with  the  pressure  transducer 
mounted  in  the  end  of  the  driven  tube,  were  concentrated  on 
investigating  the  concept  of  tailoring  to  maximize  the 
magnitude  and  duration  of  steady  pressure  in  the  reflected 
region.  While  more  than  75  shots  were  made  in  this 
configuration,  the  replacement  of  the  pressure  transducer 
and  postprocessing  of  data  taken  with  the  SD380  greatly 
enhanced  the  results.  For  this  reason,  only  data  taken  after 
run  50  using  the  SD380  are  presented.  Studies  with  the 
nozzle  attached  at  the  end  of  the  driven  tube  were 
exploratory  attempts  to  obtain  flow  visualization  in  a  shock 
tunnel  configuration.  The  results  of  the  two  series  of 
experiments  are  discussed  in  the  following  sections. 

A.  TAILORING 

Using  the  nomenclature  of  Glass  [Ref.  18],  a  tailored 
condition  is  defined  as  the  condition  when  a  Mach  wave  is 
reflected  off  the  contact  surface  and  causes  no  change  in 
the  reflected  pressure  (P5)  .  This  can  be  identified 
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experimentally  (Figure  10)  as  the  condition  when  the  ratio 
of  the  initial  reflected  pressure  (P5)  and  the  pressure 
after  the  wave  interaction  at  the  contact  surface  reflects 
off  the  endplate  (Pg)  being  equal  to  one.  Figure  31  shows 
examples  of  results  in  the  present  study  for  the  end  wall 
pressure  untailored  and  tailored  conditions  using  He  with  an 
incident  shock  strength,  M,  =  3.4.  Similar  traces  for  near 
to  a  tailored  mixture  (.7  He  and  .3  N2)  with  M,  =  2.0  are 
shown  in  Figure  32.  These  results  are  in  complete  agreement 
with  Pennelegion  [Ref.  10]. 

B.  TIME  OF  STEADY  PRESSURE 

It  was  experimentally  determined  in  the  present  study 
that  a  tailored  mixture  as  a  driver  gas  produced  a  longer 
duration  of  steady  pressure  than  that  of  a  pure  He  driver. 
The  time  of  steady  pressure  for  the  existing  NPS  shock  tube 
using  He  as  a  driver  was  approximately  3.6  milliseconds 
while  the  mixture  extended  the  test  time  duration  to  5 
milliseconds.  A  comparisons  of  the  test  times  obtained  in 
the  present  study  with  those  of  Pennelegion  [Ref.  10]  is 
shown  in  Figure  33.  The  trends  are  seen  to  be  similar. 

Since  the  work  of  Pennelegion  [Ref.  10]  was  performed  at  a 
constant  driver  pressure  of  1000  psi  in  a  N.P.L.  2"  shock 
tube  with  a  10  ft  driver  and  a  12  foot  driven  section  some 
differences  in  the  magnitudes  of  the  test  times  are  to  be 
expected. 
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Figure  31. 


HE  DRIVEa/NITRDGEK  DRIVEN 


RUN  58  IMACH  =  3.8) 


2 _  _ _ 

0  at  .OlBin  SEC  .010000 

(a)  Untailored 


Experimentally  Determined  Endwall  Pressure 
Traces  Using  Helium  as  a  Driver  Gas 
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REFLECTED  PRESSURE  RESPONSE 


7  HE  DRIVER/NITROGEN  DRIVEN  i 

RUN  70  (MACH  «  2.25) 

AT  .039ni  SEC  .020000 


(a)  Untailored 


(b)  Tailored 


Experimentally  Determined  Endwall  Pressure 
Traces  Using  a  Mixture  of  70%  Helium  and 
30%  Nitrogen  as  a  Driver  Gas 


TIME  OF  STEADY  PRESSURE  (MSEC) 


Figure  33.  Comparison  of  Steady  Pressure  Time  Achieved 
in  the  NPS  3  Inch  Shock  Tube  with  Experi¬ 
mental  Results  of  Pennelegion  [Ref.  10] 
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Fcr  low  strength  shocks  the  head  of  the  expansion  wave 
causes  an  early  degradation  of  the  reflected  pressure  in  the 
NFS  tube.  The  longer  driver  section  in  the  N.P.L.  tube  is 
responsible  for  the  difference. 

C.  PRESSURE  RECOVERY 

Figure  34  shows  the  fraction  of  initial  driver  pressure 
available  in  the  reflected  region  as  a  function  of  shock 
Mach  number.  The  data  show  that  as  shock  strength  increases 
pressure  recovery  decreases  for  both  He  and  mixed  driver 
gases.  The  data  are  in  good  agreement  with  the  results  of 
Pennelegion  [Ref.  10].  The  measured  values  of  the  reflected 
shock  pressure,  which  are  the  pressures  available  to  drive  a 
hypersonic  tunnel,  are  given  in  Tables  3  and  4. 

D.  REFLECTED  TEMPERATURE 

The  reflected  temperature  was  not  determined 
experimentally.  Instead,  shown  in  Table  2  are  values  of 
stagnation  temperatures  calculated  for  shock  strengths  in 
the  range  from  2.0  to  3.4.  These  values  were  taken  from  the 
tabulated  solutions  (based  on  equilibrium  chemistry)  of 
Bernstein  [Ref.  19:pp.  32-43]. 

E.  SHOCK  ATTENUATION 

A  representative  illustration  of  measured  shock 
attenuation  for  both  a  pure  gas  driver  and  the  driver  gas 
mixture  is  shown  in  Figure  35.  Velocity  can  be  inferred 
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Figure  34. 
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Comparison  of  Reflected  Pressure  Measurements 
in  the  NPS  3  Inch  Shoclc  Tube  with  Experimen¬ 
tal  Results  of  Pennelegion  fRef.  10] 


TABLE  3 


EXPERIMENTALLY  DETERMINED  REFLECTED  PRESSURE  MAGNITUDE 
AND  DURATION  AND  CALCULATED  TEMPERATURE  (HE  DRIVER/N2  DRIVEN) 


Run 

P5 

Time  of  Ste 

# 

Ms 

P./Pi 

(psi) 

Pe/Ps 

P5  (msec) 

50 

3.08 

994 

150 

1.79 

1.3 

51 

2.92 

353 

235 

1.23 

1.2 

52 

4 . 39 

414 

207 

1.38 

1.9 

53 

3 .32 

93 

312 

1.42 

1.5 

57 

4 . 52 

500 

183 

1.25 

0.3 

58 

3.79 

207 

304 

1.23 

0.3 

59 

3 . 66 

148 

386 

1.13 

1.6 

60 

3 . 36 

169 

264 

1.13 

3 . 6 

61 

3.72 

167 

317 

1.23 

0.5 

64 

3 .42 

164 

— 

65 

3 . 65 

158 

— 

66 

184 

67 

3.77 

173 

68 

3.81 

400 

74 

4 . 50 

667 

75 

4 . 10* 

1240 

*  Estimated  Value 


TABLE  4 


EXPERIMENTALLY 

DETERMINED 

PRESSURE 

MAGNITUDE  AND  DURATION 

AND 

CALCULATED 

TEMPERATURE 

(70%  HE, 

30%  N2 

DRIVER/N2  DRIVEN) 

Run 

P5 

Time  of  Steady 

# 

Ms 

P./Pi 

(psi) 

Pe/P^ 

P5  (msec) 

62 

2 .40 

50 

317 

1 . 10 

2 . 8 

63 

2 . 66 

51 

310 

1 .23 

1 . 5 

69 

2 .47 

50 

310 

1 . 09 

2 , 5 

70 

2.25 

50 

305 

1.11 

3 . 3 

71 

2.03 

24 

427 

1 . 02 

4 . 9 

72 

2 . 03 

22 

422 

1 . 08 

5 

73 

1.99 

21 

410 

1 . 13 

2 

77 

78 

155 

79 

3 . 0* 

171 

80 

3.4* 

638 

*  Estimated  Value 


from  the  slope  of  the  lines  between  the  thin  films  A  and  B 
and  between  B  and  the  pressure  transducer  at  C.  An 
increasing  slope  indicates  a  slower  speed.  The  general 
trend  is  that  for  a  given  gas  mixture,  the  stronger  the 
shock  strength  the  smaller  the  attenuation  between  B  and  C. 
Reflected  velocity  is  of  the  order  of  one  half  the  incident 
velocity.  Tabulated  results  for  all  runs  are  found  in 
Appendix  D. 

F.  OPTICAL  MEASUREMENTS 

Before  attempting  to  obtain  photographs  during  the 
firing  of  the  shock  tube,  a  small  freejet  from  a  1/4" 
flattened  copper  tube,  powered  by  shop  air,  was  installed  in 
the  model  location,  the  spark  light  was  adjusted  using  a 
neutral  filter  and  a  successful  shadowgraph  was  taken. 

1 .  Shadowgraph  Using  Hvcam  Camera 

The  Hycam  Movie  camera  was  triggered  manually  either 
one  or  two  seconds  before  the  diaphragm  was  ruptured.  The 
400  ft  rolls  of  ASA  250  speed  film  provided  approximately 
one  second  for  the  camera  to  accelerate,  and  a  three  second 
window  to  capture  the  event.  It  was  found  that  all  the 
available  light  was  needed  to  obtain  a  correct  exposure  and 
the  camera  was  focused  on  the  outside  surface  of  the  window. 
The  shadowgraph  was  made  visible  by  the  faint  film  which  by 
then  had  deposited  on  the  inside  of  the  window.  The  result 
was  successful. 
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with  the  blunt  body  shown  in  Figure  21(a)  installed, 
the  nozzle  starting  Mach  wave  (Figure  36(a)  and  (b) ) ,  and 
the  establishment  of  steady  flow  with  an  oblique  shock 
visible  (Figure  36(c)),  were  successfully  recorded.  Steady 
flow  (started)  conditions  lasted  a  brief  time  until  dump 
chamber  back  pressure  unstarted  the  nozzle.  Incident  shock 
strength  was  estimated  to  be  M,  =  3.8  providing  a  peak 
pressure  in  the  reflected  region  of  approximately  460  psi. 
The  actual  shock  speed  was  not  measured  during  the  test 
since  the  starting  of  the  camera  caused  premature  triggering 
of  the  spectrum  analyzer.  However  the  initial  conditions, 
and  therefore  the  reflected  pressure,  were  the  same  as  in 
run  58.  A  complete  record  of  this  run  given  in  Appendix  D. 

2 .  Spark  Shadowgraph 

Before  firing  the  shock  tube,  the  correct  exposure 
for  the  expected  shadowgraph  was  verified  by  opening  the 
film  and  firing  the  spark  source.  The  result  is  shown  in 
Figure  37.  It  was  subsequently  found  that  when  the  shock 
tube  was  fired,  the  thin  film  signal  failed  to  trigger  the 
spark,  and  thus  any  exposure  on  the  film  was  the  result  of 
light  produced  from  within  the  tube.  The  results  of 
exposing  the  film  without  discharging  the  spark  as  the  shock 
Mach  number  was  progressively  increased,  are  shown  in 
sections  (a)  to  (e)  of  Figure  38.  A  faint  glow  appeared 
first  at  a  Ms  =  3.6.  As  Mj  was  increased  for  both  a  pure  He 
driver  and  mixed  driver  gases,  the  luminosity  increased. 
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(a)  Nozzle  Starting  Shock  Wave 


(b)  Nozzle  Starting  Shock  Wave 

Oblique  Shock  V»ave 


(c)  Oblique  Shock 

Figure  36.  Results  Obtained  Using  Hycatn  Camera 
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Figure  38.  (CONTINUED) 
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Figure  38.  (CONTINUED) 


The  source  of  the  light  requires  further  investigation,  but 
it  is  clear  that  some  contamination  must  be  present  since 
the  driver  and  driven  gases,  if  pure,  would  not  be  excited  at 
the  temperature  produced  in  the  firing  process. 


65 


VI.  SIZING  CONSIDERATIONS  FOR  A  HYPERSONIC  TUNNEL 

The  experimental  program  showed  that  the  basic  shock 
tube  gave  reliable  and  repeatable  results  and,  following  the 
upgrading  recommendations  detailed  in  Appendix  A,  it  could 
well  serve  as  a  driver  for  a  hypersonic  tunnel.  Having 
experimentally  defined  the  stagnation  temperature  and 
pressure  available  to  drive  such  a  facility,  it  only 
remained  to  arrive  at  a  nozzle  and  a  dump  chamber 
configuration  suitably  matched  to  the  tube.  Simplified 
models  were  therefore  developed  to,  first,  calculate  the 
effect  of  nozzle  throat  diameter  on  the  available  steady 
flow  time  and,  second,  to  determine  a  suitable  dump  chamber 
size.  The  two  calculations  are  detailed  in  Appendix  B  and 
Appendix  C  respectively. 


A.  NOZZLE  THROAT  SIZE 

The  nozzle  throat  diameter  must  be  sized  such  that  the 
time  for  which  conditions  in  the  reflected  region  are  steady 
must  be  greater  than  the  time  required  for  the  compressed 
mass  of  gas  to  exhaust  through  the  nozzle.  The  model 
described  in  Appendix  B  estimates  the  time  in  relation  to 
the  falling  pressure  in  the  reflected  region  to  be  given  by 


t  = 


% 


C^A*P5 


1 


(B.IO) 
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If  it  is  assumed  the  supply  pressure  can  drop  to  half  its 
initial  value,  the  time  required  can  be  calculated  using 
Equation  (B.IO)  for  different  driver  gases  and  shock 
strengths.  Through  the  ratio  mo/mo,  the  calculated  time  is 
implicitly  a  function  of  the  throat  diameter.  Results  of 
such  calculations  are  shown  in  Figure  39(a)  and  (b) .  The 
four  curves  shown  illustrate  the  dependence  of  useful  time 
on  the  shock  strength  and  throat  diameter.  Clearly,  the 
shock  strength  is  less  important  than  the  throat  area. 

Based  on  these  results,  a  throat  diameter  6.5mm  (.025 
inches)  was  selected.  This  choice  made  best  use  of  the 
steady  flow  time  which  was  experimentally  determined  to  be 
from  3.5  to  5  milliseconds  using  tailored  conditions.  For 
the  6.5mm  throat.  Equation  (B.IO)  was  used  to  calculate  the 
time  for  the  reflected  shock  region  pressu_-e  to  decay  to 
various  fractions  of  its  initial  magnitude.  The  results  are 
shown  in  Figure  40.  Clearly,  for  even  the  strongest  shock, 
the  time  for  the  pressure  to  decay  to  80%  of  its  initial 
value  is  longer  than  its  experimentally  determined  tailored 
test  time  of  5  msecs. 

B.  DUMP  CHAMBER  SIZE 

The  dump  chamber  requires  a  capacity  such  that  the  back 
pressure  remains  below  the  nozzle  exit  pressure  for  the 
duration  of  steady  conditions  in  the  driven  tube.  The 
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(a)  Helium  into  Nitrogen 


Figure  39.  Time  for  Reflected  Region  Pressure  to  Drop 
to  Half  its  Initial  Value 
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TO  R 


(b)  Nitrogen  into  Nitrogen 


Figure  39  (CONTINUED) 
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analysis  in  Appenix  C  gave  the  required  ratio  of  dump 
chamber  volume  to  reflected  shock  reservoir  volume  as 


^5 


Examination  reveals  that  the  parameters  which  most  influence 
the  required  dump  chamber  volume  are  the  incident  shock 
strength  (which  determines  the  reflected  reservoir 
conditions) ,  the  desired  test  section  Mach  number,  the 
driver  gas  selected,  and  the  specified  reflected  region 
pressure  (  Pt/Pto)  at  which  the  flow  breaks  down.  Since  the 
pressure  recovery  is  a  function  of  which  driver  gas  is 
selected,  calculations  were  made  for  both  nitrogen  and 
helium.  The  results  are  shown  in  Figure  41.  A  mixture  of 
the  driver  gases  would  be  expected  to  fall  between  these  two 
curves.  Consistent  with  the  finding  above,  the  criterion 
that  the  dump  chamber  pressure  becomes  equal  to  the  nozzle 
exit  pressure  when  20%  of  the  stagnation  pressure  was 
exhausted,  was  taken  to  determine  the  necessary  dump  chamber 
volume.  Shown  in  the  figure  are  results  obtained  assuming 
the  unknov/n  temperature  in  the  dump  chamber  to  have  its 
extreme  maximum  and  minimum  values;  namely,  initial 
reflected  shock  temperature  and  initial  dump  chamber 
temperature.  From  an  examination  of  the  results,  a 


Figure  41.  Dump  Chamber  Volume  Requirements 


conservative  selection  would  be  15  cubic  meters.  This  is  a 
vessel  approximately  6  feet  in  diameter  and  15  feet  long. 

C.  HIGH  ENTHALPY  MODIFICATION 

Higher  enthalpy  flows  can  be  generated  using  the  free- 
piston  arrangement  shown  in  Figure  1.  Stalker  [Refs.  20,21] 
has  carried  out  extensive  experimental  work  in  this  area.  A 
preliminary  analysis  shows  that  the  NPS  shock  tube  has  the 
same  internal  diameter  as  Stalker's  1st  generation  free- 
piston  tunnel,  but  has  a  thinner  wall  by  1/2".  Structural 
limitations  would  prevent  the  full  duplication  of  Stalker's 
tunnel  but,  with  the  published  information  on  the  design  of 
free  piston  drivers,  it  should  be  possible  to  extend  the 
present  tube  to  an  increased  enthalpy  level  at  which  real 
gas  effects  of  interest  can  be  studied. 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  present  study  showed  that: 

The  existing  shock-tube  is  straightforward  and 
relatively  inexpensive  to  operate,  and  gives 
repeatable  test  conditions; 

The  existing  shock  tube  can  be  used  to  drive  a 
hypersonic  shock  tunnel  to  a  Mach  number  of  12 ; 

Tailored  operation  can  provide  uniform  stagnation 
conditions  for  periods  of  approximately  3.4  msec 
using  Helium,  and  approximately  5  msec  using  a  .7 
He/.3N2  mixture,  as  the  driver  gas; 

It  is  feasible  to  obtain  both  high  speed  movies  and 
short  duration  spark  flow  visualization  records  of 
the  test  section  flow  field  using  currently  available 
equipment ; 

Contamination  in  the  tube  can  lead  to  unwanted 
radiation  (which  must  clearly  be  eliminated  if  real 
gas  effects  are  the  purpose  of  the  experiment) . 

It  is  recommended  that: 

The  existing  tube  be  serviced  and  its  operating 
instrumentation  be  upgraded,  as  outlined  in  Appendix  A; 

A  Mach  10  nozzle,  15  cubic  meters  dump  chamber  and 
model  support  be  designed  and  procured  for  operation 
with  the  current  tube; 

Consideration  be  given  to  the  addition  of  a  free- 
piston  tunnel  to  provide  high  enthalpy  simulation. 

Operation  of  a  shock  tunnel  with  the  existing  tube  would 

provide  a  convenient  and  inexpensive  facility  for  use  in 

laboratory  instruction.  It  would  also  provide  a  tool  for 

thesis  studies  aimed  at  hypersonic  CFD  code  validation. 
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Extension  to  a  free-piston  driver  would  provide  an  advanced 
tool  for  hypersonic  research. 
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APPENDIX  A 


APPRAISAL  OF  THE  PRESENT  SHOCK  TUBE 


The  NPS  shock  tube  has  been  in  operation  since  the  late 
1960 's.  This  appendix  provides  a  description  of  recommended 
improvements/enhancements  of  the  existing  system,  which 
result  from  the  experience  gained  in  the  reported  program  of 
tests. 

A.  0-RINGS 

The  planned  relocation  of  the  facility  to  the  new  Gas 
Dynamics  Laboratory  requires  complete  disassembly  of  the 
shock  tube.  Many  0-rings  are  old  and  worn  and  they  should 
be  replaced  at  that  time.  Table  5  details  the  0-ring 
requirements . 

B.  THIN  FILM  GAUGES 

The  unavailability  of  replacement  thin  film  gauges  used 
to  measure  shock  speed  is  unacceptable.  Only  the  two  in¬ 
stalled  gauges  remain  operational.  The  options  are  either 
to  manufacture  thin  films  in-house  as  detailed  in  Penaranda 
[Ref.  4]  or  to  procure  them  commercially  through  the  supply 
system. 
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TABLE  5 


ORING  REPLACEMENT  GUIDE 


Location 

Number 

Required 

Width  of 
Groove 

Depth  of 
Groove 

Inside 

Diameter 

Outside 

Diameter 

(inches) 

(inches) 

(inches) 

(inches) 

Driver  Tiobe 

1 

.100 

.106 

.157 

.163 

3.502+.005 

Spaced 

Diaphragm 

2 

.100 

.106 

.157 

.163 

4.000 

4.005 

2 

.100 

.106 

.157 

.163 

3.500 

3.505 

Valve  Body 

1 

.101 

.105 

.157 

.163 

3.2485 

3.2495 

Tubing 

Sections 

3 

.101 

.105 

.157 

.105 

3.2485 

3.2495 

Instrumentation  5 
Section 

.101 

.105 

.157 
.16  3 

3.2485 

3.2495 

Driver  End  Cap  1 

.101 

.105 

.157 

.163 

3.2470 

3.2460 
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C.  PRESSURE  TRANSDUCERS 

1 .  Piezoelectric  Transducer 

The  two  Endevco  Model  2501  transducers  with  Model 
2642  M35  charge  amplifiers  should  be  returned  to  the  manu¬ 
facturer  for  calibration.  Biased  and  unbiased  responses 
should  be  checked.  Unknown,  but  apparently  unrealistic 
response  characteristics  made  these  transducers  unusable  for 
the  reported  experiments. 

2 .  Piezoresistive  Pressure  Transducer 

Procurement  of  a  piezoresistive  pressure  measurement 
system,  including  a  transducer  and  amplifier  such  as  the 
Kistler  Model  603B  Series  and  Model  5004  Charge  Amplifier, 
would  provided  continuous  measurements  of  shock  tube 
pressures . 

A  second  transducer,  to  be  used  as  a  pitot  probe  for 
flow  survey  work  in  shock  tunnel  applications,  is  also 
required.  Selection  would  be  based  on  the  selected 
nozzle  Mach  number  and  probe  geometry.  Cornell  Aeronautical 
Laboratory  publications  by  MacArther  should  be  consulted. 

D.  THIN  FILM  SIGNAL  AMPLIFIER 

The  current  thin  film  signal  amplifier  has  no  record  of 
calibration  for  the  past  ten  years.  Voltage  fluctuations 
are  apparent  and  dial  adjustment  response  is  very  poor. 

Prior  to  the  next  use  the  amplifier  components  should  be 
tested  and  replaced  as  required. 
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E.  DATA  ACQUISITION 

It  is  recommended  that  the  Scientific  Alanta  signal 
analyzer  or  its  equivalent  be  used  in  all  future  testing. 
Post  processing  capability,  and  the  availability  of  four 
channels  with  independent  trigger  adjustment  provide  a 
valuable  experimental  capability  not  available  with  the 
Tektronics  549  and  551  oscilloscopes. 

F .  DIAPHRAGMS 

The  existing  brass  diaphragms  allow  driver  operation  at 
a  constant  600  psi  level.  Redesign  to  a  higher  operational 
level  of  approximately  700  psi  would  increase  the  range  of 
reflected  pressures  available  to  drive  a  shock  tunnel,  with 
an  associated  cost  of  higher  gas  consumption.  While  not 
essential,  this  modification  would  reduce  problems  which 
would  be  expected  in  an  extension  to  a  shock  tunnel  with 
respect  to  the  measurement  of  small  aerodynamic  forces. 

Gaydon  [Ref.  2:p.  91]  and  Glass  [Ref.  18:  pp.  503-511] 
provide  valuable  design  insight  into  successful  diaphragm 
material  selection  and  preparation. 

G.  MISCELLANEOUS  HARDWARE 

Valve  7  and  valve  11  located  in  the  driven  end  of  the 
tube  [Ref.  4 : Figure  2]  exhibited  some  signs  of  leakage 
during  the  loading  process.  Repair  or  replacement  of  the 
valve  seats  is  recommended. 


79 


APPENDIX  B 


ESTIMATION  OF  SHOCK  TUNNEL  NOZZLE  THROAT  AREA 


In  order  to  size  a  suitable  nozzle  to  be  driven  by  the 
high  pressure  and  temperature  in  the  reflected  region,  a 
simplified  model  was  developed.  The  model  was  based  on  the 
following  assumptions: 

A  constant  volume  reservoir  corresponding  to  the  in¬ 
itial  reflected  shock  conditions  was  available  to  drive 
the  nozzle  flow.  The  reservoir  volume  was  defined  by 
the  diameter  of  the  tube  and  the  axial  distance  from 
the  endplate  of  the  tube  to  the  point  where  the  contact 
surface  and  reflected  shock  met  (Fig.  3) ; 

The  gas  remaining  in  the  reflected  region  expands 
isentropically  to  fill  the  constant  volume; 

The  flow  through  the  nozzle  throat  was  isentropic  and 
always  choked.  (This  implies  that  the  evacuated  dump 
chamber  had  sufficient  volume,  so  as  not  to  influence 
the  nozzle  flow. ) 

The  total  mass  initially  in  the  reflected  reservoir  is 
given  by. 


where 


^0  R  T 


V 


g  5 


0  ' 


(B.i: 


^2 
rD  X 


^^0  = 


The  mass  at  any  later  time  is  given  by 


m 


R  T. 

g  t 


V, 


(B.2) 


As  the  mass  flows  through  the  choked  nozzle  the  rate  of 
change  of  mass  in  the  volume  is  equal  to  that  which  flows 
out  through  the  nozzle,  thus 


dm 

dt 


A*P, 
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^  g 

2 
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Integrating 


+-  t  A  *P 

-  r  dm  =  f  — ^  c,  dt  (B.4) 

0  0'  vT^ 

Since  the  expansion  in  the  reservior  is  assumed  to  be  isent- 
ropic,  the  stagnation  pressure  and  temperature  are  related 
according  to 


T 


5 


(B.  5) 


Taking  the  ratio  of  Equations  B.l  and  B.2  and  substituting 
into  Equation  B.5, 
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Substitution  into  Equation  B.3  gives 


(B.6) 


dm 

dt 


so  that 
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On  integrating,  the  variation  of  mass  in  the  reflected 
reservior  with  time  is  given  by 


t 


/ _ 2_ 

CiA*P5  \^Y5-1 


(B.9) 


and  the  variation  of  pressure  in  the  reflected  region  with 
time  is  given  by 
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To  estimate  a  suitable  nozzle  throat,  Equation  B.IO  was 
programmed  and  a  ratio  of  Pt/Pj  equal  to  one  half  was 
selected  as  the  criterion  for  calculating  the  useful  run 
time.  Using  throat  area  as  a  variable  all  other  parameters 
were  a  function  only  of  the  incident  shock  strength  and  of 
the  physical  size  of  the  tube.  Thus  Equation  B.IO  enabled 
the  available  time  for  steady  flow  through  a  given  throat 
diameter  to  be  estimated. 
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APPENDIX  C 


ESTIMATION  OF  THE  REQUIRED  DUMP  CHAMBER  VOLUME 

In  a  shock  tunnel  application,  the  dump  chamber  pressure 
is  required  to  remain  below  the  nozzle  exit  static  pressure 
for  the  useful  duration  of  the  ref lected-shock  supply  condi¬ 
tions  in  the  shock-tube.  To  aid  in  sizing  a  dump  chamber 
for  a  Naval  Postgraduate  School  shock  tunnel  the  lollowing 
analysis,  based  on  a  simplified  model,  was  developed.  The 
model  assumes  the  pressure  and  temperature  are  uniform  in 
the  dump  chamber,  and  that  the  velocity  there  is  effectively 
zero.  It  is  assumed  that  the  flow  to  the  exit  of  the  nozzle 
is  isentropic,  and  that  the  volume  of  the  nozzle  is  negligi¬ 
ble  in  comparison  to  that  of  the  chamber.  The  volume  of  the 
dump  chamber  is  given  in  terms  of  the  initial  mass  contained 
and  its  properties,  by 
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RTi^itio 


D 
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Dividing  by  the  reflected  reservoir  volume,  similarly  writ¬ 
ten  in  terms  of  the  initial  mass  contained  and  its  proper¬ 
ties,  the  ratio  of  dump  chamber  to  reflected  shock  volume  is 
given  by 
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As  the  flow  continues,  the  original  reflected  pressure  (P5) 
decreases  to  lower  values  (Pt)  .  Since  the  nozzle  area  ratio 
is  fixed,  the  pressure  at  the  exit  of  the  nozzle  (P^)  also 
steadily  decreases.  Meanwhile,  the  pressure  in  the  dump 
chamber  (P^)  steadily  increases  until,  at  a  point  in  time, 

Pd  =  Pg,  and  the  useful  test  time  is  over.  The  ratio  of  the 
initial  reflected  pressure  to  dump  chamber  pressure  is  then 
given  by 
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A  mass  balance  gives 
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so  that 
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Using  Equation  B.6  from  Appendix  B 
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and  substituting  in  Equation  C.6, 
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Using  Equation  C.2  applied  to  the  initial  conditions. 
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Combining  Equations  C.2  through  C.9 
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Equation  C.IO  provides  a  relationship  between  Vd/Vq  and  the 
value  of  Pt/p5  which  terminates  the  useful  test  time,  in 
terms  of  given  initial  reflected  reservoir  and  dump  chamber 
conditions.  Using  an  interval  halving  technique  the  dump 
chamber  volume  requirements  were  investigated  using  an 
assumption  that  the  temperature  in  the  dump  chamber  was 
equal  to  its  initial  value,  and  also  with  the  worst-case 
assumption  that  it  was  equal  to  Tj.  (It  is  noted  that  the 
temperature  change  can  be  included  explicitly  through  the 
conservation  of  energy.  However  the  complexity  of  doing  so 
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is  not  warranted  when  an  estimate  of  the  needed  volume  is 
all  that  is  required.) 
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APPENDIX  D 


SUMMARY  OF  EXPERIMENTAL  DATA 

This  appendix  contains  data  taken  using  thin  film  gauges 
and  the  pressure  transducer.  First,  in  Figure  42,  the 
measurements  obtained  of  shock  Mach  number  using  oscillo¬ 
scopes  (proceeding  Run  50)  are  shown.  For  all  runs  from  Run 
50  on,  data  were  taken  using  the  Scientific  Atlanta  SD380 
signal  analyzer.  Tables  6  and  7  list  data  for  runs  using  a 
mixture  of  helium  and  nitrogen  and  for  runs  using  pure 
helium  as  a  driver,  respectively.  For  each  run,  trace  A 
recorded  the  response  of  the  thin  film  at  location  A,  trace 
B  recorded  the  response  of  the  thin  film  response  at  station 
B,  trace  C  recorded  the  response  of  the  endplate  pressure 
transducer.  Delay  times  were  programmed  prior  to  each  run. 
Traces  B  and  C  were  delayed  by  the  same  time.  Figure  43 
provides  definitions  of  items  which  appear  on  the  SD380 
screen  and  identifies  the  menu  selections  of  interest. 

The  last  35  pages  of  the  appendix  contain  printouts  of 
the  SD380  screen  for  the  runs  listed  in  Tables  6  and  7. 

Each  printout  can  be  interpreted  with  reference  to  Figure 
43.  Times  were  obtained  for  shock  travel  between  trans¬ 
ducers  on  traces  A,  B  and  C,  using  the  time  scale  given  on 
the  printout,  and  the  delay  times  listed  in  Tables  6  and  7. 
These  results  are  included  in  the  two  tables. 
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Figure  42.  shock  Strength  vs.  p,/p,  for  Data  Reduced 
from  Polaroid 


TABLE  6 


TIME  MEASUREMENT 

(70%  HELIUM  30%  NITROGEN  DRIVER/N2  DRIVEN) 


Tempera¬ 

Delay 

Delay 

Time  A 

Time  B 

Time  C 

Run 

ture 

A 

B 

to  B 

to  C 

to  B 

# 

(F) 

(msec) 

(msec) 

(msec) 

(msec) 

(msec) 

62 

63 

-1.01 

0.00 

2.84 

1.67 

3 . 63 

63 

65 

-1.01 

-1.01 

1.77 

1.77 

1.72 

69 

64 

-1.01 

0.49 

1.81 

— 

70 

64 

-1.01 

0.51 

1.36 

1.83 

3.86 

71 

64 

-1.01 

0.51 

1.72 

2 . 03 

3 . 79 

72 

65 

-1.01 

0.51 

1.76 

2 . 03 

3.82 

73 

64 

-1.01 

0.51 

1.80 

2.07 

3.86 

TABLE 

7 

TIME  MEASUREMENT 

(HE 

DRIVER/N2 

DRIVEN) 

Tempera¬ 

Delay 

Delay 

Time  A 

Time  B 

Time  C 

Run 

ture 

A 

B 

to  B 

to  c 

to  B 

# 

(F) 

(msec) 

(msec) 

(msec) 

(msec) 

(msec) 

50 

66 

0.00 

0.49 

1.45 

51 

66 

0.00 

0.49 

1.64 

0.94 

2.50 

52 

62 

-1.01 

0.49 

0.92 

1.24 

2.48 

53 

62 

-1.01 

0.49 

1.30 

0.92 

2 . 55 

57 

62 

-1.01 

0.49 

0.92 

1.09 

2 . 50 

58 

62 

-1.01 

0.49 

1.19 

1.13 

2.48 

59 

62 

-1.01 

0.49 

1.17 

1.23 

2 . 34 

60 

62 

-1.01 

0.49 

1.32 

1.11 

2.44 

61 

62 

— 

— 

1.15 

— 

64 

64 

-1.01 

0.49 

1.23 

— 

65 

63 

-1.01 

0.49 

1.30 

66 

61 

— 

... 

67* 

63 

— 

— 

1.18 

68* 

63 

— 

— 

1.18 
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0.51 

0.97 
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Figure  43.  SD380  Menu  Selection  Items  of  Interest 
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APPENDIX  E 


COMPUTER  PROGRAM  FOR  SHOCK  TUNNEL  EVALUATION 

The  attached  program  was  used  for  the  evaluation  of  the 
NPS  shock  tube  and  shock  tunnel.  Sources  of  shock  tube 
relationships  are  noted  in  the  program  were  applicable. 
Output  was  tabulated  by  the  program  and  graphical  analysis 
was  provided  by  the  plotting  routine  'EASYPLOT'  installed  on 
the  NPS  IBM  computer. 
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FILEi  SH0CK2  FOR^RAM  A1 


CC  THIS  PROGRAM  CALCULATES  THE  IDEAL  PERFORMAIICE  OF  A  SHOCK 

CC  TUBE  USItIG  THE  ANALYSIS  PRESEtUED  IN  THE  'FUNDAMENTALS 

CC  OF  GAS  DYNAMICS'  BY  J.  A.  OWCZAREK  PG  395-398  (FIG  8-22) 

8  'THE  SHOCK  TUBE  IN  HIGH  TEMPERATURE  CHEMICAL  PHYSICS'  BY  GRA 
GRAYDON  HURLE  1963. 


IMPLICIT  REAL  *  A  (A-H,M,P-Z) 

REAL  FCN, XI , X2 . XTOL , FTOL , FCNl , FCN2, FCN3, FCNA 
INTEGER  I, NLIM,J, 11,112 
EXTERNAL  FCN, FCNl , FCN2 , FCN3 , FCNA 

DIMENSION  M1(200),M2(200),M3(200),PAP1(200),P2P1(200) 

DIMEtlSION  P3PA(200),MS(200),US(200),T2(200),P5P2(200),P5FA(200) 
DIMENSION  C3(200),CA(200),P4P3(200),T3(200),TC(200),XC(200) 
DIMENSION  VC(200),P5PI(200),T5T1(200),UR(200),DELTIM(200) 
DIMENSION  DIST (200),TCC( 200 ),XC1(200), DELTA (20 0),PDI(8) 

DIMENSION  XMS(200),XP2P1(200),XP5P2(200),XP3PA(200),XR02R1(200) 
DIMENSION  XT2T1(200),URUS(200),TSTART(200),VDVO(8),VD(8) , AA(8) 
DIMENSION  P5(200),T5(200),DTHR(200),DT(200),MT(200),DTEST(200) 
DIMENSION  PT(200) ,TT(200),T(200, A),TMA5S(200),DTDTHR(5),P1(200) 
DIMENSION  FT(A,10),PTPTO(10) 

COMMON  /B/  R1,MM,C8,DTHR,T 

COMMON  /A/  CONST, GAHMA1,GAMMA4, AAAI  SO, Cl, Cl SQ, CCS, CCA, C5, AAAI, C9 

COMMON  /C/  ATASTR,MT,DTEST 

COfIMON  /E/  DTHR,T,MT,DTHRDT,MM,C8,R1 

COMMON  /D/  P5,T5,T1,AA,PDI,J, 11,112 

COMMON  PTPTO,FT 

DATA  XI, X2, XTOL, FTOL, I, NLIM/0. 1,6. 0,0. 00 01, 0.0 0001, 0,50/ 


SELECT  GAS  CONSTANTS 
GAS  MW  GAMMA 

AIR  28.97  l.A 

ARGON  39.94  1.67 

HE  A  1.67 

N2  28.02  l.A 

C02  AA.Ol  1.29 


SPEED  OF  SOUND 
345  M/S 
321 
1062 
550 
269 


DRIVEN  GAS  CONSTANTS 

Rl=831A/28.02 

Tl=288 

GAMMA1=1 .A 

Pl=l . OEOA 

DRIVER  GAS  CONSTANTS 

RA  =  851A/28 .02 

TA=288 

GAMMAA-1 . 40 

PA=600X6 .89A8E03 

A  IS  INCREMENT  ON  PAPl  LOOP 

A  =  100 

MM  IS  THE  NUMBER  OF  LOOPS 
MM  =  6 

PAP1(1)=20 


INPUT  TUNNEL  PARAMETERS! DRIVER  LENGTH  LI  8  LA) 
Ll=6 . 453 
LA=2.375 
DT1=3X . 0254 

C  TO  EVALUATE  TEST  CONDITIONS  SET  FL2=1 

THR1=.0065 
FL1  =  0 

C  FL2  CALCULATES  TEST  SECTION  MACH  f 

FL2  =  1 

C  FL3  IS  FOR  NOZZLE  SIZING  CALCULATIONS 

FL3  =  0 

C  CALCULATES  TEST  SECTION  CONDITIONS 

FLA  =  1 
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C  FL6  EXPERIMENT  OUTPUT  RUNS 

FL6  =  0 

C  FL7=I  PRINTS  MISC  RELATIONSHIPS 

FL7  =  0 

C  FL8=1  EVALUATES  PT/PTO  RATIOS  FOR  GIVEN  THROAT 

FL8  =  1 

C  EVALUATE  SPEED  OF  SOUND  CONSTANTS 

A1=SQRT(GAMMA1*R1*T1) 

A<4  =  SQRT(GAMMAAXRA»T4) 

A4A1=A4/AI 
A4A1SQ=A4A1XX2 


EVALUATE  CONSTANTS  USED  COMMONLY  USED 
C 1  =  ( GAMMA !  +  !)/( 2  XGAMMAl ) 

ClSQ  =  Cl3f*2 

C5=(-2!(GAMMA4)/(GAMMA4-I) 

C6  =  ( GAHMA4-1 )/ ( 2*GAMMA4 ) 

C7  =  ( GAMMA4  +  1 )/ ( 4*GAMMA4 ) 

C8=(GAMMA1-1 )/(2*GAMMAl)*( 1 ) 

C9  =  (GAMMA1  +  1)/(2)«(GAMMA1-I)  ) 
C10=(GAMMA1/(GAMMA1-1) ) 


THIS  LOOP  INCREMENTS  PRESSURE  P4/P1  BY  A  AND  THEN 
USES  INTERVAL  HALVING  TO  DETERMINE  M5 . 

DO  1  J=1,MM,1 
C0MST  =  P4P1(  J) 

CALL  INTHV(FCN,X1,X2,XR,XT0L,FT0L,NLIM,I) 

M3( J)=XR 
C3( J)=CC3 
C4(J)=CC4 
X1=XR 
X2=6.0 
1  =  0 

THE  ROOT  HAS  BEEN  FOUND  AT  XR  IF  I  EQUALS  1  OR  2 
K  =  J  +  1 

P4P1(K)=F4P1( J)+A 
I  CONTINUE 

THIS  LOOP  USES  P4/P1  8  M3  TO  DETERMINE  P2 , PI , MS , M2 , US , P3P4 

DO  5  1=1, MM, I 

FIND  P2/P1  (EQN  8-57) 

P2P1(I)=P4P1(I)/C4(I) 

M2(I)=S0RT(29f(P2Pl(I)  +  (l/P2Pl(n)-2)/(GAMMAlx(GAMMAl-l)* 
&((GAMMA1+1)/(GAMMA1-1)+P2P1(I)))) 

FIND  MS  a  SPEED  OF  SHOCK  (E0N(8-51) 

MS(n=SQRT((P2Pl(I)-l)XCl  +  l) 

US(I)=MS(I)XA1 
FIND  P3/P4 
P5P4(I)=C3(I)X)fC5 
P4P3(I)=1/P3P4(I) 

EVALUATE  MAXIMUM  DURATION  OF  UNIFORM  FLOH  IN  THE  REGION 
BETWEEN  THE  SHOCK  WAVE  AND  THE  CONTACT  SURFACE. 

REFEREtICE  PAGE  435  PROBLEM 

T2(I)  =  {2/(P5P4(I)X*C7)-(4*A4)(Cl-r3P4(I)K5fC6)*(l/P3P4(I))!<xC7)/ 
8((GAMMA4-l)*Al»S<3RT(l  +  (GAMMAl  +  l))f(P2PlCI)-l)/ 
8(2*GAMMA1)))))(L4/A4 


T3  IS  THE  THEORITICAL  MAXIMUM  DURATION  OF  UNIFORM  FLOW 
IN  THE  REGION  BETWEEN  THE  CONTACT  SURFACE  AND  THE  TAIL  OF 
CENTERED  EXPANSION  WAVE  (REGION  5).  (PG  433  PBLM  8-6) 

T3(1)=(L4/A4)X((GAMMA4-1)/2)*((P4P3(I)xx(3-GAMMA4)) 

8*X(1/(GAMMA4X4)))/(1-(P3P4(I)XKC6)) 


TC  a  XC  ARE  THE  TIME  AND  DISTANCE  AT  WHICH  THE  WAVE  FORMED  BY 
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THE  REFLECTION  OF  THE  CENTERED  EXPANSION  HAVE  OVERTAKES  THE 
CONTACT  SURFACE.  (REF  PGA33  PBLM  8-7) 

TC(I)  =  (L<i/AA)x2!‘-(PAP3(I)K*C7) 

XC(I)=(L4X4/(GAMMA4-1)K(1-(P3PA(I)KXC6)/(PAP3(I)KKC7))) 


CALCULATION  OF  REFLECTED  REGION  PARAMETERS  ACCORDING  TO  REF  2 

VELOCITY  OF  CONTACT  SURFACE  VC 

VC(I)=2)(A1/(GAMMA1  +  1  )»(MS(I)-1/MS(I )  ) 

PRESSURE  AND  TEMPERATURE  IN  REFLECTED  REGION  5  ( EQN  2.31 
P5Pl(I)  =  ((2>fGAMMAl);MS(I)**2-(GArirtAl-l  ))/(GAI!MAl  +  l  )  ) 
&*((  (3*GAMMA1-1)«MS(  I)»«»2-2»<(GAMMA1-1)  )/( (GAMMAl-l)?fri3CI  )*>!2  +  2)  ) 
T5T1(  I)  =  ((23«(GAMMAl-l  ) XMSt  I ) X?; 2+ (  3-GAMf1Al )  ) 
83«((3XGAMMA1-1  )xMS(I)«5«2-2X(GAMMA1-1  ) )/ ( ( GAtlMAl  +  l)  2*113  ( I  )  ?<  V2  )  ) 
CALCULATE  VELOCITY  OF  REFLECTED  SHOCK  UR  (EQN  2.38) 
UR(I)=US(I)»({2+(2/(GAMMAl-l))!((l/P2PlCI)))/ 
&((GAMMA1+1)/(GAMMA1-1)-1/P2P1(I))) 

OBSERVATION  TIME  AND  DISTANCE  INREFLECTFD  SHOCK 

DELTIM(I)  =  L1K(US(I )-VC(I))/(US(I )*(URCI)+VC(  I) ) ) 
DIST(I)=DELT1M(I)kUR(I) 

TIME  TO  EXPANSION  WAVE  REACH  CONTACT  SURFACE  (EQN  4.7) 
C=((GAMMA4  +  1)/(23«(GAMI1A4-1)  )  ) 

TCC<  I)=2*L4*(1  +  (GAMMA4-1  ) *M3 ( I  )/'2 )  ^^C/AA 
XCK 1 )=TCC(I)K((2XA1*(MS(I )K*2-1 ) )/( (GAMMAl  +  l )*MS(  I ) ) ) 
TEST  TIME  BETWEEN  CONTACT  SURFACE  AND  SHOCK 
DELTA(I)=  TCC(I)*(1-2J((MS(I)XK2-1 )/( (GAMMAl+1 )*M3( I )**2) ) 
P5P2(I)=P5P1(I)/P2P1(I) 

5  CONTINUE 


SECOND  METHOD  OF  CALCULATIONS 
Xl=.l 
X2  =  6  .  0 
XTOL=0.0001 
FT0L=0. 00001 
1  =  0 

NLIM=50 
XR=  5 

DO  2  J=1,MM,1 
C0NST=P4P1( J) 

CALL  1NTHV(FCN1,X1,X2,XR,XT0L, FTOLrNLIM, I) 

Xri3(  J)=XR 

XP2F1( J)=(2*GAMMA1KXMS( J)»*2-(GAMMA1-1) )/(GAMMAl+l ) 

XP3P4( J )  =  1/(P4P1( J)/XP2P1( J  ) ) 

XR02R1( J)=(<GAI1MA1+1)KXM5( J)**2)/((GAMMA1-1 )*XM5( J)**2>2) 
XT2T1( J )=( (GAnnAl*XMS( J)*X2-(GANriAl-l )/2 ) ) X( ( GAMMAl -1 )/2*XM3( J )** 
SE+D/C  ((GAMMAl  +  1  )/2  )*X2*XMS(  J  )x.*2  ) 

XP2Pl(J)  =  <2*GAHMAl*nS(J) *5^2- (GAMMA  1-1 ) )/(GAnnAl  +  l ) 

XP5P2( J)  =  (( (GAMMA1  +  1)/(GAMMA1-1) )  +  2-l/xr2Pl( J) )/ 
8(1+(GAMMA1+1 )/(GAMMAl-l)*l/XP2Pl( J)) 

XP3P4( J)=1/(P4P1( J)/XP2P1(J)) 

XR02R1( J )=( (GAMMAl+1 )*MS( J )*K2 )/ ( ( GAMMAl -1 ) »MS( J ) **2+ 2 ) 
XT2T1( J )=( (GAMMA 1*MS(J)*X2-(GAMMA1-1 )/2 ) ) *( ( GAMMAl-1 )/2«M5( J ) ** 
82+l)/( ( (GAMMA1+1)/2)x*2*MS(J)**2) 

X1=XR 
X2=6  .  0 
1  =  0 

M3  HAS  BEEN  DETERMINED  BY  INTERVAL  HALVING 
DATA  X1,X2,XT0L,FT0L,I,NLIM/XR, 5. 0,0. 0001, 0.00001, 0,50/ 

URU5( J)=UR( J)/US( J) 

P5P4( J)=P5P1( J)/P4P1(J) 

C  P5( J)=P5P1( J)XP1 

F5( J)=P5P4( J)»P4 
Pl( J)=P4/P4P1(J) 

T5( J)=T5T1( J)XT1 

WRITE(*,»)  P5P1(J),P4P1(J),P1( J),UR( J),DELTIM( J) 

2  CONTINUE 

C  129 


FILE:  SH0CK2  FORTRAN  A1 


C=====^-r====^^==r=-=======================::=====r =============== 

IF  (FL3.EO.O)GO  TO  36 

C  THIS  SECTION  DETERMINES  P5  T5  AND  DETERMINES  TESTIME  AS  A 

C  FUNCTION  OF  NOZZLE  THROAT  DIAMETER. 

NN  =  I 

DO  3  1=1, MM 

C  REFLECTED  MO  AVAILABLE  TO  DRIVE  NOZZLE 

C  =  SQRT(GAMMAl/Rl)*(l-fGAMMAl-l  )/2 ) ( -1  >!C9  ) 
TMASS(I)=P5(I)»(UR(I)*DELT1M(I)>X3. 1A17*DT1XX2/(R1»(T5(I))(A) 

C  HRITE(x,x)  P5P1(I),P5(I).T5(I),NN,TMASS(I) 

CALL  N0ZZLE(P5(I).T5(I),NN, TMASSC I) ,C, FL3) 

C601  FORMAT! 'PAP1=*,D9. 3, 2X, •P1=*,D9. 3, 2X, 'TOTAL  MASS  MO= ' , D9 . 4, 2X, 

C  a'RUN  =',15/) 

NN=NN+1 

3  CONTINUE 

DO  309  1=1,100 

HRITE(8,602)  DTHR< I ) , T( I , 1 ) , T ( I , 2 ) , T ( I , 3 ) , T ( 1 , 4 ) 

602  F0RMAT(5(E9.3,2X)) 

C201  FORMAT! 'P5=' , E9 . 2 , 2X, ' T5= • , E9 . 2, 2X ,  'NN='  , II 
C  «,2X, 'TOTAL  MASS  M0=',E9.2) 

309  CONTINUE 
C  WRITE!8,603T 

C603  FORMAT!/, 'THROAT  DIAM' ,2X, 'RUNl ' ,6X, 'RUN2' ,7X, 'RUN3' ,7X, 'RUN4' ,/) 
C  IIRITE!8,601)  ! P4PI ! I) , PI ! I ) , TMASS ! I ) , I , I = 1 , 4 ) 

C  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

C  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

36  IF  !FL8.EQ.0)G0  TO  37 

C  THIS  SECTION  DETERMINES  P5  T5  AND  DETERMINES  TESTIME  AS  A 

C  FUNCTION  OF  PT/PTO  FOR  A  GIVEN  THROAT  DIAMETER 

NN  =  1 

DO  82  1=1, MM 

C  REFLECTED  MO  AVAILABLE  TO  DRIVE  NOZZLE 

C  =  SQRT!GAMMA1/R1)X!1-!GAHMA1-1 )/2  )  X* ! -1 XC9  ) 
TMASS!I)=P5!n*!UR!l)xDELTIM!  I  )  )X5. 1417XDT1XX2/!R1XT5!  I  )X4) 
CALL  N0ZLE2!P5!I),T5!I),NN,TMASS!I),C, FL3) 

HN=t(N+l 
82  CONTINUE 

DO  89  1=1, MM 

IIRITE!8,802)  PTPTO! I) ,  FT!  1 , 1 ) ,  FT !  1 , 2) ,  FT  !  1 , 3 ) ,  FT !  1 , 4 ) 

802  F0RMAT!5!E9.3,2X)) 

89  CONTIflUE 

HRITE!8,803)  ! P4PI ! I ) , 1=1 , MM) 

803  F0RMAT!/,4! 'P4P1=',E9.3,2X)) 

C=  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  =  = 

37  IF  !FL2.Eq.O)GO  TO  133 

C  DETERMINE  TEST  SECTION  MACH  »  AS  FUNCTION  OF  THROAT 

C 

C  CALL  EXPAN!THR1,DT1) 

C  SUBROUTIflE  EXPAN!THR1,DT1) 

C  REAL  X1,X2,XT0L,FT0L,DT1,THR1 

C  INTEGER  I,NLIM,J 

C  DIMENSION  MT!200).DTEST!200),DTDTHR!5) 

CC  COMMON  /A/  CONST, GAMMA1,GAMMA4,A4A1S0, Cl, C150,CC3,CC4,C5,A4A1 

C  COMMON  /C/  ATASTR,MT,DTEST,C9,DTDTHR 

HRITE!8,300) 

DTEST!1)=DT1 
Xl=.10 
XR  =  50 

DO  30  J=l,8 
X2=100 
XTOL  =  0 . 0001 
FTOL=0 . 00001 
1=0 

NLIM=50 

ATASTR=!DTEST! J)/THR1)»X2 
C  WRITE!*, X)  THR1,DTE5T!J),ATASTR,C9 

CALL  INTHV!FCN3,X1,X2,XR,XT0L,FT0L,NLIN, I) 

MT!J)=XR 
X1  =  XR 
X2=100 
C 

K=J  +  1  130 
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500 

301 

30 


DTEST(K)=(DTEST(J)+DT1X.5) 

DTDTHRC J)=DTEST(J)/PT1 
PRINT  'HELP' 

HRITEC8,301)  MT( J ) , DTDTHRI J ) , DTEST ( J ) , ATASTR, THRl , DTI 
FORMAT (/, AX, • MTEST • , IX, ' DTEST/DTUBE' , IX, * DTEST ' , 2X, ' AT ASTR ' 
8,AX, 'DTHROAT* ,3X, 'DTUBE'  ) 

F0RMAT(A(F9.3),2(2X,E9.3)) 

CONTINUE 


133  1F(FLA.E0.0)G0  TO  105 

THIS  SECTION  EVALUATES  TEST  CONDTIONS  FOR  VARIOUS 
TEST  CHAMBER  AREAS  GIVEN  A  FIXED  NOZZLE  THROAT  DIAMETER 
OUTER  LOOP  IS  FOR  PAPl  CHANGE,  INtJER  LOOP  EVALUATE  CH  AREA 
DO  22  J=1.MM 
IIRITECS,  400) 

DO  21  1=1,8 
HRITE(K,»)  CIO 

C  IIRITECX.AOS)  MT(I).P5CJ),T5(J),PT(I),TT(I),C10 

C4  06  FORMAT (2X,  'MACHSeX,  'PS  =  ' , 6X,  •  PT  ' ,  5X ,  •  T5  '  ,  8X,  '  TT  '  ,  5X,  'TSTART'  ) 

C405  FORMAT(6(D9.2)) 

CALL  TCOND(MT(I),P5(J),T5(J),PT(I),Tr(I),TSTART(I),THRl,C10) 
IJRITE(8,4  01)  MT(I),P5(J),PT(I),T5(  J),TT(I),M5(  J),P4P1(  J),PH  J) 
8,TSTART(I) 

400  FORMAT!/, 'TEST  MACH ’, 5X, ' P5 *, 4X, ' TEST  P ' , 4X , • T5 ' , 4X , 

8 'TEST  TEMP' ,3X, 'MS' ,4X, ' P4P1 • . 6X, ' PI ' , 7X, ' 1  ST  ART ' ) 

401  F0RMAT(9(E8.3,1X),E6.1) 

21  CONTINUE 

C  DUMP  CHAMBER  SIZING  EVALUATION 

l/RITE(8,411) 

DO  461  112=1,1 
DO  460  11=1,1 
C  AA(I1)=.1XI1 

AA(I1)=.8 
PDI<I1)=. 013332 
Xl=-10000 
XR=50 
X2=100000 
XTOL=  .  0001 
FTOL=. 000001 
1  =  0 

MLIM=200 

C  WRITECS,*)  AA(I1),PDI(I1) 

CALL  IMTHV(FCN4,X1,X2,XR,XT0L, FT0L,NL1M, I ) 

VDVOCIl )=XR 

C  MRITE(8,X)  AA(I1),PDI(I1),VDV0(11) 

X1=XR 
X2  =  10 

VD( II )=VDVO( II )»(UR( J)XDELTIM( J )»3. 1417XDT1XX2/4) 

WRITE! 8, 412)  P5! J ) , T5( J ) , VDVO! 1 1 ) , VD! 1 1  ) , AA ! 11 ) 

411  FORMAT !/,3X, ' P5 * , 8X, ' T5 ' , 8X, • VDVO ' , 7 X , ' VD ' , 8X , 'PTPTO' , 

86X) 

412"-  F0RMAT!5!E9.2,1X)) 

460  CONTIflUE 

461  CONTINUE 

HRITE!8,402)  DELTIM! J ) , T! 1 , J  ) 

402  FORMAT!/, 'TEST  TIME  OF  SHOCK  TUBE  =  ',D9.4,/ 

8, 'TIME  FOR  .5XP5  OUT  OF  NOZZLE  =',D9.4) 

C  8, 'TIME  TO  START  NOZZLE  lEHPERICAL  EST IMATE ) ' , D9 . 4 ) 

22  CONTINUE 

105  K=MM 

IF  !FL6.EQ.O)GO  TO  699 
HRITE!8,695) 

WRITE! 8, 696)  ! P4P1 !J ) ,MS{ J ) , US! J ) , VC ! J ) , URUS ! J  ) , T5 ! J  ) 
8,P5!J), J=1,K) 

695  FORMAT! 2X, ’XP4P1 ' , 7X, ’MS' , 9X, 'VS' , 9X, 'VC , 9X, 'URUS', 

85X, 'P5' ,8X, ’T5'  ) 

696  FORMAT!6!F9.3,2X),E9.3) 


699  IF  !FL1 .EQ.O)GO  TO  108 

WRITE!8,201)  131 
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HRITE(8,20  0)  ( PAPl ( J ) , MS( J ) , XF2P1 ( J ) , XP5P2 ( J ) , XP3P4 ( J  ) 
8.XR02P1( J),XT2T1( J),XMSCJ),URUS(J^ , J=1.K) 

200  F0RMAT(9(F11.3,2X)) 

201  FORMAT (2X, 'XPAPl •.7X, 'MS ' , 8X, • XP2P1 ' ,8X, • XP5F2 ’ , 6X, • XP3PA ' 
a,‘iX,  •R0U2/R0U1',4X. •XT2/T1',6X, •XMS',6X, *  UR/US' ) 

IIRITE(8,203)  (P4P1( J),MS(J),XT2T1( J),T5T1(J),P2P1( J) 
&,P5P^(J),DELTIM(J),J=1,K) 

203  FORMAT(7(D9 .4,2X)) 

108  WRITE(8,204) 

IIRITE(8,205)  R1,T1,GAMMA1,R<4,T4,GAMMA4 

204  FORMAT! 3X, ’P4P1 ',8X, 'MS',8X, 'XT2T1 ' , 7X,  ' T5T1 * , 7X,  ’P2P1  • 
&,7X, 'P5P4',5X, 'DELTIHE') 

205  FORMAT! 'Rl=', F9. 2, 2X, 'Tl^* , F9 .2, 2X, 'GAMMA1= ' , F9 . 2, 

8/, 'R4=',F9.2,2X, 'T4=*,F9.2,2X, 'GAMMAA^ • , F9 . 2) 

C=======r==3=============================:=== 

IF  !FL7.EQ.O)  GOTO  109 
C  OUTPUT  DATA  TO  A  FILE 

11RITE!8,103) 

IIRITE!8, 100)  !P4P1!I),HS!I),M2!I),M3!I),A4A1, 

8PSP1!I ),T5T1!I),P3P4!I),P2P1!I),P5P2! I),DELT1MII), I=1,K) 
IIRITE!8,106) 

WRITE! 8, 107)  ! P4P1 ! I ) , UR! I ) . DIST ! I ) , XCl ! I ) , TCC! I ) , DEL T IMI I ) , 
8DELTA!I),T2!I),1=1.K) 

10«  F0RMAT!4X, 'P4P1 • ,2X, ’VEL  REFL ' , 4X , ' DI S  REFL',4X, 

8'EXP  P0S'.4X,'EXP  TIMES4X. 'DEL  TIM ' ,  SX ,  '  DEL  T  A  '  ,  8X ,  '  T  2  •  ) 
WRITETS, 104  ) 

107  F0RMAT!8!F9.4,2X)) 

WRITE! 8, 105)  ! P4P1 ! I ) , T2! I ) , T3! I ) , TCI  I ) , XCC I ) , 1  =  1 , K ) 

100  F0RMAT!11!F9.4,1X)) 

102  FORMAT !4X, ' P4P 1 ' , 8X , • T2 ' , 8X, »P3P4',6X, 'P2P1 ' ) 

103  FORMAT !4X, ’ P4P1 • , 6 X , 'MG  * , 

89X, 'N2',6X,  M13',8X,'A4A1',6X,'P5P1' 

&  6X  *  T5T 1  *  6X 

six, •P3r4' ,6X, 'P2P1' ,6X, •P5P2',6X, • TESTTIME ' ) 

104  FORMAT  !2X,  'P4P1  M2X,  ^2*  ,8X,  •T3M0X,  '  TC  * ,  9X ,  '  XC '  ) 

105  F0RMAT!5!F10.4,2X)) 

109  X=1 
STOP 
END 


REAL  FUNCTION  FCNIX) 

REAL  X 

COMMON/ A/  CONST, GAHMA1,GAMMA4,A4 Also, Cl , Cl SQ, CCS, CC4,C5,A4A1 ,C9 
CC3=! 1+IGAMMA4-1 )»X/2) 

CC4  =  !CCS*X!2*GAMMA4)  X^-Il/IGAMMAA-l)  ) 

FCN!=!  1  +  !  .  5!<A4A1S0)X!  !GAHHA1XX/CC3)X5(2) 

8XIC1  +  SQRT! !ClSQ  +  4/!A4AlSQX!GAMMAlXX/CCS)X)f2)  )  )  )  ) 

8XCC4-C0NST 

WRITE! 8, 200)  FCN,GAMMA1,GAMMA4, A4A1 50, CCS, CC4, CONST 
200  FORMAT!7<F8 .2)  ) 

RETURN 

END 

REAL  FUNCTION  FCNIIY) 

REAL  Y 

COMMON/ A/  C0NST,GAMMA1,GAMMA4, A4A1SQ,C1,C1S0,CC3,CC4,C5, A4A1,C9 
FCN1=! !2XGAMMA1*YX)«2-GAMMA1+1)/!GAMMA1+1 ) )*! l-I IGAMMA4-1 )/! A4A1 
8XIGAMMA1+1 ) )X!Y-1/Y) ))XKC5-C0NST 
PRINT?!,FCN1,C0NST,A4A1,C5 
RETURN 
END 


REAL  FUNCTION  FCH3IZ) 

REAL  Z 

COMMON  /C/  ATASTR,MT,DTEST 

COMMON  /A/  CONST, GAMMA1,GAHMA4.A4A1S0, Cl, CISQ, CCS, CC4, C5, A4A1 ,C9 
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C  HRITECK,*)  ATASTR,GAMHA1,Z,C9 

FC!^3  =  ATA?TR-1'Z*<2“(1-^(CAMMA1-  1  )/2»Z*w2) 

8/((GAMMAl  +  l))))(*C9 
C  HRITE(x,*)  ATASTR,GAMMA1,Z,FCN3 

RETURN 
END 
C 

REAL  FUNCTION  FCNAd-J) 

REAL  II 

C  INTEGER  J,I2. 11 

DIMENSION  P5(200),t5(200),AA(8),PDI(8),MT(200) 
C0MI10N/C/  ATASTR,MT.DTEST 
COMMON/ D/  P5,T5,T1,AA.PD1, J.I1,I12 

FCN9=N-(l+PDiai)yT5(J)»IJ/(P5{J)>!Tl  ) -AA ( 1 1  ? .  7 1  4 ) ^Tl  » 
8(1+.2*100)*X3.5/(AA(I1)XT5(J)) 

WRIT E< 8, 131)  AA(I1),PDI(I1),P5( J),H,T5(J), J,I1 
131  FORMAT ( 'A= ' , E9 . 2 , • PDI = • , E9 . 2 , * P5= ' , E9 . 2 , 
8'lt=',E9.2,lX,E9.2,lX,F4.1,lX,F4.1) 

RETURN 
END 


SUBROUTINE  NOZZLEC X. Y, NN,MT0TAL , C0N5T2 , FL 3 ) 

THIS  ROUTINE  DETERMINES  TIME  FOR  NOZZLE  FLOH  FOR  A  GIVEN  THROAT 
DIAMETER. 

REAL  )(8  X,Y,MT0TAL,C0NST2,P1 
INTEGER  NN,FL3 

DIMEMSIOtl  DTHR(20  0),T(200,4) 

COMMON  /A/  C0NST,GAMMA1,GAMMA4, A4A1SQ,C1 ,C1SQ,CC3,CC‘4,C5,  A4A1,C9 

COMMON  /B/  P5,T5,R1,MM,C8,DTHR,DTHRDT,T 

COMMON  /B/  R1,MM,C8,DTHR,T 

FI=3. 1416 

DTHR(1)=.0065 

DO  30  1=1,100,1 

T(I,NN)=MT0TALXDSQRT(Y)X2X(1/.5xkC8-1 )/(C0NST2xri kDTHRC I ) **2/4 «X 
8*(GAMMA1-1)) 

J  =  I+1 

DTHR( J)=DTHR(I)+.0005 
IIRITE(8,200)  DTHR(I),T(I,NN) 

200  F0RMAT(2(E9 .3,2X) ) 

201  FORMAT?  'P5=',E9. 2, 2X,  ’T5=’ ,E9 .2,2X,  'NN=M1 
S,2X, 'TOTAL  MASS  M0=’,E9.2) 

30  CONTINUE 

IIRITE(8,201)  X,Y,NN,MTOTAL 

RETURN 

END 


SUBROUTINE  NOZL E2( X, Y, NN , MTOT AL , C0NST2 , FL 3 ) 

THIS  ROUTINE  DETERMINES  TIME  FOR  NOZZLE  FLOW  FOR  A  GIVEN  THROAT 
DIAMETER  HITH  VARATION  OF  PT/PTO 
REAL  *8  X,Y,MTOTAL,CONST2,PI 
INTEGER  NN,FL3 

DIMENSION  DTHRC200),T(200,4),PTPTO(10),FT(4,10) 

COMMON  /A/  C0NST,GAMMA1,GAMMA4, A4A1SQ,CI,C1SQ,CC3,CC4,C5,A4A1,C9 

COMMON  /B/  P5,T5.R1,MM,C8,DTHR,DTHRDT,T 

COMMON  /B/  R1,MM,C8,DTHR,T 

COMMON  /F/  PTPTO,FT 

PI=3.1416 

DTHR(1)=.0065 

PTPT0(1)=.2 

DO  80  1=1,4 

FTC I,NN)=MT0TAL*DSQRT(Y)*2X(1/PTPT0(I)**C8-1 ) 

8/(C0NST2*PI*DTHR( 1 )**2/4XX 
a*CGAMMAl-l )) 

J  =  I  +  1 

PTPT0(J)=PTPT0(I)+.2 
80  CONTINUE 
RETURN 
END 

SUBROUTINE  EXPANC THRl , DTI ) 
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C  REAL  X1,X2,XT0L.FT0L,DT1,THR1 

C  I1.7EGER  I,NLIM,J 

C  DIMENSION  MT(200),DTEST(200),DTDTHR(5) 

C  COMMON  /A/  CONST, GAMMA1,GAMMAA,AAA1SQ, Cl, CISQ,CC3, CCA, C5,A<4A1,C9 

C  COMMON  /C/  ATASTR,MT,DTEST,C9, DTDTHR 

C  DTEST(I)=DT1 

C  Xl=I.O 

C  XR  =  2 

C  DO  30  J=l,2 

C  X2=10 

C  XTOL=0.001 

C  FTOL=0.0001 

CC  1  =  0 

C  NLIM=50 

C  ATASTR=(DTEST( J)/THRl)x»2 

C  HRITEC*,*)  THR1,DTEST{J),ATASTR 

C  CALL  INTHV(FCN3,X1,X2,XR,XT0L,FT0L,NLIM,I) 

C  MT(I)=XR 

C  X1=XR 

C  K=J+1 

C  DTEST(K)=(DTEST(J)+DT1) 

C  DTDTHR( J)=DTEST( J)/DT1 

C  WRITE{8,301)  MT( J),DTDTHR(J),DTEST( J) 

C300  F0RMAT(2X, * MT EST ’ , 2X , • DT EST/ DTUBE * , 2X , 'DTEST' ) 

C  MRITE(8,300) 

C301  F0P.MAT(3(F9.2)) 

C  30  CONTINUE 

C  RETURN 

C  END 

C 
C 
C 

SUBROUTINE  TCONDC Z1 , P5, T5, PT, TT , TIME , THRl , Cl 0 ) 

REAL  X8  P5,T5,TT, TIME, Z1,THR1,B, CIO 

COMMON  /A/  CONST, GAMMA1,GAMMAA,A4A1SQ, Cl ,C1SQ,CC3  'C4 , C5, AAAI , C9 
C  COMMON  /£/  DTRH(200),T(200),Mf,DTHRDT,MM,C8 

COMMON  /B/  Rl,MM,C8,DrHR,r 
HRITECX, A03) 

HRITE(X,A0A)  Z1,P5,PT,T5,TT,C10 
MRITE(X,X)  CIO 
A5=S0RT(GAMMA1*R1XT5) 

PT=P5/(l+(GAnMAl-l )/2»ZlXX2)KXC10 
TT=T5/(1+(GAMMA1-1 )/2*Zl*K2) 

C  B  IS  THE  ANGLE  OUT  OF  THE  THROAT  30  DEGREES 

B=TAN( 15*3. 1A16/180) 

C  TIME=0 

TIME=.76*THR1/(2*B*A5)*((D5QRT(Z1) )*(1-1/Zn 
8+l/15*((DS0RT(Zl) )**3-l)  +  l/25*((DS'3RT(Zl)  )**5-l)) 

A03  FORMAT! 2X, 'MACH' ,6X, 'P5  = ' , 6 X, ’ PT ' , 5X , ’ T5 ’ , 8X , ' TT ' , 5X, ' T STAR T ' ) 
MRITECX, AOA)  Z1,P5,PT,T5,TT, CIO, TIME 
AOA  FORMAT! 7 (D9. 2)) 

RETURN 

END 

C  - 

C 

C 

c 

c 

c 

SUBROUTINE  INTHV! FCN, XI , X2, XR, XTOL , FTOL , NL IM, I ) 

C 

C  - 

C  THIS  ROUTINE  MAS  DEVELOPED  BY  PROF.  MILLER  DEPT  AERONAUTICS 

C  SUBROUTINE  INTHV  : 

C  THIS  SUBROUTINE  FINDS  THE  ROOTS  OF  A  FUNCTION  , 

C  F!X)  =  0,  BY  INTERVAL  HALVING. 

C 

C  - 

c 

C  PARAMETERS  ARE  i 

C  FCN  -FUNCTION  THAT  COMPUTES  VALUES  FOR  FIX). MUST  BE  DECLARED 

C  EXTERNAL  IN  CALLING  PROGRAM. IT  HAS  ONE  ARGUMENT,  X. 
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FILE:  SH0CK2  FORTRAN  A1 


XI,  X? 
XR 

XTOL, 

NLIM 

I 

1  =  1 
1  =  2 
I=- 
I=- 
HHEN 
PRINT 


-INITIALS  VALUES  OF  X.F(X)  MUST  CHANGE  SIGN  AT  THESE  PTS 
-RETURNS  THE  ROOT  TO  THE  MAIN  PROGRAM. 

-TOLERANCE  VALUES  FOR  X,F(X)  TO  TERMINATE  ITERATIONS. 
-LIMIT  TO  NUMBER  OF  ITERATIONS. 

-A  SIGNAL  FOR  HOW  ROUTINE  TERMINATED. 

MEETS  TOLERANCE  FOR  X  VALUES. 

MEETS  TOLERANCE  FORF(X) 

NLIM  EXCEEDED. 

F(XI)  NOT  OPPOSITE  IN  SIGN  TO  F(X2) 

THE  SUBROUTINE  IS  CALLED, THE  VALUE  OF  I  INDICATES  WHETHER  TO 
EACH  VALUE  OR  NOT. 1=0  MEANS  PRINT  THEM,I.NE.O  MEANS  DON'T, 


FTOL 


REAL  FCN, XI, X2,XR, XTOL, FTOL 
INTEGER  NLIM,I,J 
REAL  F1,F2,FR,XERR 


CHECK  THAT  F(X1)  8  F(X2)  DIFFER  IN  SIGN. 

FI  =  FCN(Xl) 

F2=  FCN(X2) 

IF  (F15(F2  .GT.  0.0)  THEN 
I  =  -2 

PRItIT  201 
RETURN 
END  IF 


COMPUTE  SEQUENCE  OF  POINTS  CONVERGING  TO  THE  ROOT 

DO  20  J=1,NLIM 
XR  =  (Xl+X2)/2.0 
FR  =  FCN(XR) 

XERR  =  ABS(Xl-X2)/2.0 
IF  (I  .EQ.  0.0)  THEN 
PRINT  199,  J,XR,FR 
END  IF 


CHECK  ON  STOPPING  CRITERIA 

IF  <  XERR  .LE.  XTOL)  THEN 
1  =  1 

PRIfIT  202,  J,XR,FR 
RETURN 
END  IF 


IF  (  ABSCFR)  .LE.  FTOL  )  THEN 
1  =  2 

PRINT  203,  J,XR,FR 
RETURN 
END  IF 


COMPUTE  NEXT  X 


20 

C 


IF  (  FR!<F1  .GT.  0.0  )  THEN 
X1=XR 
F1  =  FR 
ELSE 
X2  =  XR 
F2  =  FR 
END  IF 
CONTINUE 
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HHEN  LOOP  IS  NORMALLY  COMPLETED, NLIM  IS  EXCEEDED 
I  =  -l 

PRINT  200,  NLIM,XR,FR 
RETURN 


199  FORMATC  AT  ITERATION* , 13, 3X, '  X= * , El 2 . 5 , AX , '  F(X)  =',E12.5) 

200  FORMAT!/'  TOLERANCE  NOT  MET  AFTER  ',IA,'  ITERATIONS  X  =', 

1  E12.5,'  F(X)  =  •,E12.5) 

201  FORMAT!/'  FUNCTION  HAS  THE  SAME  SIGN  AT  X1&X2') 

202  FORMAT!/'  TOLERANCE  MET  IN  ',IA,'  ITERATIONS  X  =',E12.5, 

1  '  F!X)  =  ',E12.5) 

203  FORMAT!/'  F  TOLERANCE  MET  IN  ',IA,'  ITERATIONS  X  =  ',E12.5, 

1  '  F!X)  =',E12.5) 

END 
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